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ON TRIPLETS OF PLANE CURVILINEAR 
ELEMENTS WITH A COMMON SINGULAR POINT 


By C. C. HSIUNG (Chekiang) 
[Received 21 April 1946] 


1. Introduction 
Ir is known that, if two plane curves have a point O for their 


common singular point of the same kind with the same tangent, so 
that the power-series expansions for non-homogeneous projective 
coordinates representing these two curves in the neighbourhood of 


a y=ark+... and y= bak+... a #5), 

then the ratio a/b is a projective invariant of the two curves. 
C. Segre* gave a simple geometrical characterization of this invariant 
for any value of k. About four decades later E. Bompianif inter- 
preted it again, by using certain algebraic curves, for the case in 
which & is a fraction. Since no corresponding projective invariant 
can be determined for two curves having a common singular point 
of different kinds, it is desirable to undertake the investigation for 
three curves. This is the object of the present paper. 


2. Derivation of an invariant 

Suppose that C,, C,, C, are three plane curves having O as a 
common singular point of different kinds but with the same tangent, 
so that the behaviours of the common tangent ¢ at the point O with 
them are distinct. Let the homogeneous projective coordinates of a 
point in the plane be denoted by (2, y, z). If we choose the point O 
to be the vertex (0, 0, 1) and the common tangent ¢ to be the side 
y = 0 of the triangle of reference, then the power-series expansions 
of the curves in the neighbourhood of O can be written in the forms 


(1) 
(2) 


(3) 


* C. Segre, ‘Su alcuni punti singolari delle curve algebriche, e sulla linea 
parabolica di una superficie’, Rendi. dei Lincet (5), 6 (1897), 168-75. 

+ E. Bompiani, ‘Alcuni risultati di geometria proiettivo-differenziale’, 
Rendi. del Semin. Mat. e Fis. di Milano, 10 (1936), 1-28. 

3695.18 K 
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Without loss of generality we may assume the exponents A, yp, v in 
(1), (2), (3) to be distinct integers or fractions greater than 1. 

In order to find a projective invariant associated with the singular 
point O of the curves C,, C,, C,, I consider the most general projective 
transformation of coordinates that leaves invariant the point O and 
the common tangent ?: 

L = Ay, 2'+A,py’, J = Gey’, Z = Ag, X'+Agy'+Gg32'. (4) 
The effect of this transformation on equations (1), (2), (3) is to pro- 
duce three other equations of the same form whose coefficients, 


indicated by accents, are given by the formulae 


a.,a\510’ = ada, Ay, a4, 1b’ = at, b, Beg 83 'c’ = a,c. (5) 


Further, elimination of @,;, G9, @33 from equations (5) shows imme- 
diately that the expression 

I = at—h’1-# (6) 
is a projective invariant associated with the singular point O of the 
curves C,, Cy, C5. 


3. Geometrical characterizations of the invariant J. 

In this section I shall characterize geometrically the invariant J 
obtained in the previous section. It is convenient to distinguish four 
cases, according as A, pw, v are integers or fractions. 

Case I. All three exponents X, 1, v integers. Let the vertices (1, 0, 0), 
(0, 1, 0) of the triangle of reference be denoted by O,, O, respectively. 
It is easily seen that, associated with the singular point O of the 
curve C,, we can uniquely determine an algebraic curve of order A 
having a (A—1)-ple point at O, with all tangents coinciding in the 
line O, O,, and having a contact of order A with the curve C, at O. 
From (1), the equation of this curve is readily found to be 

yz-1_ax = 0. (7) 

From (2), (3), we can likewise obtain the equations of two analogous 

algebraic curves of orders p, v associated with the singular point O 
of the curves C,, Cs, 

1.e. yz+-1— bat = 0, (8) 

y2”-1—cx” = 0. (9) 

Elimination of z from (7), (8) yields the equation of the line joining 

the point O to any intersection, other than O, of the curves (7), (8): 


qe-1 1/(A—p) 
ae (3) y= 0. (10) 
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Likewise, the line joining the point O to any intersection, other than 
O, of the curves (7), (9) is given by the equation 


ar-1 1/(A-v) 


Thus we obtain the following geometrical characterization: 
T'he cross-ratio of the line OO,, the common tangent t, and the two lines 
(10), (11) as equal to 


—1\ 1/A—p) —1\ 1/(A-v) 
(0. - ais # (5 " — JA-D)A-pyA-»), (12) 





pA-1 on 
Case II, All three exponents A, p, v fractions. For this case we 
may put 
ya p=, yor td (13) 
Az Pe Ws 
where A, > Ag, Hy > Me, Vy > Va, and (Aj, Az), (41, He), (vy, ve) are three 
pairs of relatively prime integers. First of all we observe that, asso- 
ciated with the singular point O of the curve C,, we can uniquely 
determine an algebraic curve of order A, satisfying the following 
conditions: 
(i) its polar curves of O, of orders 
1,..., Ay—Agq—1, Ay—A_+-1,..., Ay—1 
are indeterminate, while that of order A,—A, degenerates into 
the line O, O, counted A,—A, times; 
(ii) O is a singular point at which the curve has contact of order 
A,—1 with ¢; 
(iii) the curve has contact of order A, with C, at O. 
From A = A,/A, and (1), it is easy to obtain the equation of this 
algebraic curve, namely 
yr zi —aha = 0. (14) 
If the curves C,, C, are used instead of C,, two similar algebraic 
curves of orders p4;, v, can be determined whose respective equations 
are 
ye zha-H2— bee gh — 0, (15) 
y"22-"—c%r a" = 0, (16) 
Elimination of z from (14), (15) shows immediately that the equation 
of the line / joining O to any intersection, other than O, of the curves 
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(14), (15) can be written as (10). Similarly, the line /’ joining O to 
any intersection, other than O, of the curves (14), (16) is given by (11). 
Thus we arrive at the following geometrical characterization of the 
invariant J for this case: 

The cross-ratio of the line OO,, the common tangent t, and the two 
lines l, l’ above is equal to [4-DIA-wA-), 

Case III. Two only of the exponents X, p, v are integers. Let A, uw 
be integers and v a fraction denoted by v,/v, as in the previous case. 


Now the three curves (7), (8), (16) can be used for our purpose. If p 
be the line joining the point O to any intersection, other than O, of 
the curves (7), (16), then the cross-ratio of the line OO,, the common 
tangent t, and the two lines (10), p is equal to [0-VIA-wyxA-), 


Case IV. Only one of the exponents X, uw, v is an integer. Let A be 
an integer and yp, v fractions denoted by p/p», v;/vo as in Case I. 
If g, g’ be the lines joining O to two intersections, other than O, 
of the curve (7) with the curve (15) and with the curve (16) respec- 
tively, then the cross-ratio of the line OO,, the common tangent t, and 
the two lines q, q' is equal to [A-DIA-wA-»), 








ON SYSTEMS OF ALGEBRAIC EQUATIONS AND 
CERTAIN MULTIPLE EXPONENTIAL SUMS 


By 8. H. MIN (Petping) 

[Received 10 July 1946] 

1. THROUGHOUT the paper we use p to denote a prime. Let 
f(x) = a, a*+...4a,2 (k > 2), 
be a polynomial with integer coefficients. The exponential sum 
Pp + 
S{f(x)} seen > e2rif(z)[p | 
z=1 

is called a generalized Gaussian sum. Mordell proved that, if a, is not 
divisible by p, then, as p tends to 00, 

Sif(a )} = O(p'- Wk) 
where the constant implied by O depends only on k. 


The n-dimensional generalization of S{f(x)} is the n-ple exponen- 


tial sum 
p p (2 


GI e,,...,%_)} = Foe Fe expy— Se beiad Tn) }> 


%=1 t=1 
where f(2,,...,%,,) is a polynomial of degree k with integer coefficients. 
More generally we consider a field K containing p” elements and 
suppose that f(z,,...,x,) is a polynomial of degree k in the field K. 
Construct the sum 


SEF ey), K} = Lex e[feey tall), 


where G[a] is the trace of a, and 2j,...,”,, run respectively over all 
elements in K. 
Now let , 
F (0, t4-Bys-++5 pn E+Bn) = 2 Flay. ns Bys+++s Bye” 
and 
(2k, oR ok, of) 
Oy , , Oxy ‘ , 6 n | 
wae. s | ee "ea. - 
af, oh, oF, of, 
ory . Coty op, , eB, 
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I shall prove that, if k > 2n and M is of rank 2n, i.e. if a certain 
2n < 2n minor is not identically zero, then, as p tends to infinity, 


y ‘ , i _— —I/k)\ 
S{f (aty,--+)%p), K} = Ofpmna-Uny, 


where the constant implied by O depends only on k and n. 

When m = n = 1, M is of rank 2 unless f(x,) reduces to a linear 
function of x,, so that the theorem includes Mordell’s theorem as a 
special case. In fact, when n = 1, the minor determinant formed by 
the first two rows is 


Sail f'(Bs)f"(B1)—2{F"(Bs)}"I. 


This is identically zero only when f(x,) is a linear polynomial. 

When x = 2, the theorem was established by Hua and me. We 
also proved by an elaborate method* that, if M is of rank < 2n = 4, 
f(x,,%,) can be transformed by a linear non-singular transformation 
in K to g(x,y) which is either of degree < $k in y or is a polynomial 
in #?+7y? (r€ K). From this we deduced that the theorem holds 
unlessf f(%,, x2) is a polynomial in a single variable which is itself a 
linear function of x and y. It is very likely that a similar result exists 
in the general case, but this seems to be very difficult to establish. 
In the two-dimensional case I have obtained a still better result, 
which, roughly speaking, may be regarded as a two-dimensional 
analogue of a result due to Davenport.{ [hope to publish it elsewhere. 

2. In the theory of a system of linear equations the most important 
role is played by the matrix of the system. A natural generalization 
of this is the ‘Jacobian matrix’ of a system of algebraic equations. 


Let Ai se (§.== 0.590) (1) 


* ‘On a double exponential sum’, Science Report of Tsing-Hua Univ. (in 
the press). A short sketch of the proof has been published in the Science Record 
of Academia Sinica, vol. i, Nos. 1-2. It is impossible, however, to give a 
sufficiently complete short sketch. 

+ The method of procedure where f(x, 7.) is a polynomial in x?+-7y? may 
be outlined as follows (taking m = | for simplicity). The sum is the same as 
the corresponding sum with 2?+7y? replaced by x+y, provided that each 
term is multiplied by 

(1+ (/p)}{1+ (y/p)}- 
This gives four sums; only the one containing (zy/p) presents difficulty. This 
is dealt with by putting y = xz, and applying Mordell’s theorem to the resulting 
sum over x. 

t ‘On certain exponential sums’, J. fiir die reine und angew. Math., 169, 
Heft 3 (1933), 158-76. 











ON EXPONENTIAL SUMS 135 


be a system of algebraic equations in a field F. We write the 
Jacobian matrix as 







‘th, | 
OL y 

> | me ns F a 3 ‘i | (2) 
| Ox, Ox, J 


and, when m = n, denote its determinant by J(f,,...,f;n)- 


DEFINITION. A solution of (1) for which the matrix ¥(f;,..-,fm) 18 of 
rank r is called a solution of rank r. If r = n, the solution is said to be 
non-singular ; otherwise it is singular. 


The theorems to be proved hold both for fields of characteristic 
zero and those of characteristic p. We are interested in large primes p 
and polynomials of fixed degree. 

Corresponding to the fact that a system of linear equations has 
only one solution if its determinant is different from zero, is the 


following 


THEOREM 1. The number of non-singular solutions of (1) is O(1), 
where the constant implied by O depends only on n, m and the degrees 


of f, (¢ = 1,...,m). 


Proof. Without loss of generality we may assume that m = n. 
The theorem is evident for n = 1. Suppose that the theorem is true 
for n—1 variables and n—1 equations. Let us prove the theorem for 
n variables and n equations. 

If f,,.--,f, do not contain* x,, the theorem is true, for in this case 
J (f,,--,f,) = 0 identically and (1) has no non-singular solutions. If 
only one of them contains z,, the theorem is also true; for supposing 
that f,, contains x, and that (x, , 2) ,, 2) is a non-singular solution 
of (1), then we aes 


ae co ‘ee £0 teu 


which shows that (x{,..., 2 ,) is a non-singular solution of f; = 9,..., 
f,-1 = 0 and that f,( 210) +, ,,z,,) is not identically zero. Let the 
jae of f; with respect - x, be A;. Then the theorem is true when 

* We consider /,,...,f, as polynomials in x, whose coefficients are poly- 


nomials in 2,...,%,_,- But f; may be of degree 0 in x,. In this case we say, 
for convenience, that f; does not contain 2,. 
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A, +...+A, <1. For, in this case, at most one of the equations 
(1) contains x,. Now suppose that the theorem is true when 


Ait... +A, <7. 

Let us prove the theorem when d,+-...+-A, = r+. 

Without loss of generality we assume that 

fi frotet+..+huay fa = fanxr +.-+for, 

when A, > A, > 1 and f,;,f2; are polynomials in 2j,..., 2,4. 

Consider the system of equations 

fi = fof: —frot fa = 0,* f,=90 (0 =2,...,m). (3) 

Evidently, every solution of (1) is a solution of (3). I shall prove that 
each non-singular solution (2{°,...,2') of (1), for which fy) 4 0, is a 
non-singular solution of (3). 

By differentiating the first equation in (3) with respect to 2;, 


n), we have 


foi fom oF fs 4 Foe p, — Se tah, 


= faz ~ —fio arde—te c at x, = 2 (i = 1,...,m—]), 


vi 


of . af , —Ne af 2-1 
of, = fos 1 —fig2® “ —(A,—Ag) fig tH -** If, 


OX c "4 n 





= Soa —f; 9 t- nits at Ly = 
v OX», 


vy, 
Therefore 
nae . ¢ . oa ae 4 
I (fisSe-sbn) = food (fishin) FO at x, = x7”. (5) 
Now suppose that (2{,..., 2) is a non-singular solution of (1), for 
which f., = 0. Then the solution satisfies both 
= Ree As 
f= 0, fs = feo—Sa0 Xn = 0, 
and fi = 9, foo = 9, fi 
By (6), 
fo fs fro » Oe 
- —_ (. = l,...,.»—1), 
Ox; On; Ox; 
ofs = fe a yhr-l ofe at x. = 7 
ox, Ot, °* n% 0x, ‘ R 
* In the case in which this equation becomes an identity, (3) has no 
singular solutions. 
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Hence 


OAS fret) = I fobofaSd+0r I frfrorfasfn) (8) 
at x, = 2, This shows that J(f,fafgs--rf,) and I(fysfoorFar-+-fa) 


cannot both be zero. 

Therefore, by (5) and (8), each non-singular solution of (1) must 
be a non-singular solution of one of the systems (3), (6), (7). Since 
the sum of the degrees (with respect to x,,) of each of the systems is 
less than or equal to 7, the theorem follows by induction. 

I shall establish three theorems which have some interest in them- 
selves but are not required for the treatment of the exponential sums. 


THEOREM 2. The number of solutions of (1) of rank r can be arranged 
in O(1) sets,* each set being such that, if a certain n—r of the variables 
are given, the remaining variables are determined with only O(1) 
possibilities. 


Proof. For a solution of rank r, at least one of the r-rowed minors 
in 3(f,,---:fm) is not zero. Let it be the minor in the left upper corner, 
say M. By Theorem 1, if 2,,,,...,%, are given, (1) has only O(1) 
solutions for which M ~ 0. The theorem follows immediately. 

3. Now let us consider homogeneous equations 


Gl Xq»---»Fnis) =O (6 = 1I.,...,m), (1) 


where the g),...,9,, are homogeneous polynomials in 2},...,%,4;- Con- 
sider the matrices 














[ 91 | 

OX, , OXy +1 
MI: prery Im) — ? 

9m Gm 

eee ml 

(a * 3%... oe 

Ot, Ogg Bbgyy Ol n+1 
3j(Gas-++» Gm) = ° . ° ° . ° ° . ° . ° . . = 

Gm 9m Gm =: Gm | 

LOX, OX;_4 OX 544 OXy +1 


We have the following theorems. 


* The constant implied by O in this theorem and Theorem 4 depends only 
on m, n and the degrees of the equations considered. 
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THEOREM 3. There are no solutions for which 
Rank of 3; < Rank of 3 (j = 1,2,....n+1). 
except the trivial solutions 
4%, =090, ..., Bay = O. 
Proof. If k; is the degree of g; and x, = 2? (v = l,...,.n+1) isa 
solution of (1), then 
Og; Cq; 
0— kg, = a, t+ me ee. * at 2, = 2, 
OX, OX n+1 
This shows that the rank of 3 is equal to that of 3; at x, = 2{ unless 
x\° — 0. The theorem follows. 

THEOREM 4. T'he number of non-equivalent solutions (proportional 
solutions being called ‘equivalent’) of rank r can be arranged in O(1) sets, 
each set being such that, if a certain n—r of the variables are given, the 
remaining variables are determined with only O(1) possibilities. 

Proof. Consider first the solutions with x; ¢ 0. These solutions 
satisfy also 


- x Ln4 . ’ 
25 *9 (2450.2) F_43) = 02 ) =) (¢ = i...) (1’) 
x; x; 

J J 
where k; is the degree of g; (1 <i <m). Let X; = x;,/a; or x;,,/2; 


according as i < j—l ort >j. Then (1’) can be written as 








y as - ° ” 
Gh A eyry hg) =O (6 = I....,00). (1") 
Plainly 
[2G G1 | 
|oX, eX, 
. 
| avy ary 
OG» OG» 
— « © + mae 
(eX, ex, J 
(91 1 _ a G1 
| - ae - < eee Se 
| Cary OX;-1 Oz 44 OX 
— gr kitt | 
a 
| Ym Im 9m 9m 
Ox, OX» Oj 4) OL» 
— a~kjt+l 
= ag Mth ye 


As shown in the’ proof of Theorem 3, for a solution with 2; ¢ 0, the 
rank of 3; is equal to that of 3. Hence, corresponding to a solution 
(#1,-+-)%y4,) Of (1) of rank r with x; ¢ 0, there is a solution (X,,..., X,) 
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of (1”) of rank r. The theorem therefore follows from Theorem 2 if 
we consider only non-trivial solutions. But, since there is only one 
trivial solution, the proof is thus completed. 

4. Before going on to prove the result stated in §1, I establish 


two lemmas. The constant implied by O in this section depends 
only on k and n. 


LemMA 1. Let K be a finite field containing p™ elements, where p is 
a prime and m a positive integer. Let 
f(x) = a, a*+...4a, (a, 40, a; Ek). 
Then > e27iS[f(x)\/p ae O(pma-Uh), (1) 


where S{a] denotes the trace of a, and x runs over all elements of K.* 
Proof. Let @ be the generating element of K with respect to the 
ground field 7 consisting of the residue classes to modulusp. Then 
every element x of K can be written as 
X= ata O+...pa, 10"? (a; € 7). 
Since S[a+y] = S[x]+ Sly], we have, forae K, 


™ Pp Pp — = = 
> e27iSlazl/p — > 3 pa E27 i(Xo Sa] +..+%m-1 S[ad”"))/p 
ao=1 Xm-1 =1 


_ {pm if S[a*]=0 for k= 0,...,m—1, 
~ \0 otherwise. 
Further, S[a6*] = 0 for k = 0,...,m—1 impliesa = 0. In fact, we 
have Slax] = 0 for all ae K. If a 4 0, we should have 
0 = Sfaa"] = 1, 


which is absurd. Hencet 


s pas ri S[b1 04.040, 2" I!p 2k 


3 bz 


= > ~ > > ae > > net >») eR ES (Dy foe HEY 1 one Ya) tee HOM Hoe —Y DID 
Zz Te U1 


VE by be 


M &s 


Je 3 z : S pm _ O(p?"*). (2) 
71 Te V1 os Uk 
Lyte FLEA Hoe FU 
b+. -yt+..+9 


* Setting m = 1 and G[x] = x we get Mordell’s theorem. In the proof we 
can simply put 6 = 0 and omit all unnecessary words. 


+ We consider 4j,..., a, as fixed and },,..., 6, as variables that assume all 
values in K. 
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If A and p» are both in K and A = 0, then 


> e2riSfAr+piip — > e27iS[M(a)\/p 
Zz < 


The transformation x = Ax’+p (Ac K, »e K, A ~ 0) transforms 
the sum in (1) into at least p”(p”— 1)/k different terms in the sum (2). 
In fact, if f(Aw+p) = f(x),* then 

Ma, = ay, m1 ya,,+-AF-la,_ = Ay_y- 
The first equation gives at most k values of A and the second equation 
determines p uniquely. Hence 
re | > e27iS[K(x)l/p ak — O(p?”*), 


and > ertiSimip — O(pmi-ub), 
7 

LEMMA 2. Let f(x,,...,x,) be a polynomial not zero in K, then the 
number of solutions of 
| J F feys.0.5%,) = 0 (3) 
as O(p™-»). 

Proof. The proof is evident for n = 1. Suppose that the lemma 
is true for n—1 unknowns. Let f(2},...,2,) = gpa*+...1+g; where g; 
are polynomials in 2,...,x,_, and g, # 0. The number of solutions 
of f = 0, gy ~ 0 is evidently O(p™"-) and that of f = 0, g, = 0 is, 
by our supposition, O(p™"-»p”) = O(p™-»). Therefore the number 
of solutions of (3) is O(p™"-»). 

Now we can prove 

THEOREM 5. Let f(2},...,x,) be a polynomial of degree k > 2n in the 
field K containing p™ elements. Let 


k 
I ( a | b-+-B 0005 Xn t+B,,) = 2 Elo Xs By,--r Ba 


Suppose the matrix 








[2h oF, ) 
OX OP n 
a | 
ap ap | 
oF, oF, | 
Oxy eB, J 
is of rank 2n, i.e. a certain 2n x 2n minor is not identically zero. Then 
ba yee > e2riS[f(x1 —_ z,)I/p — O(pmra-Ub), 
at In 


* This means that they are the same polynomial in K. 
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Proof. Let J(a,,...,8,) be a minor which is not identically zero. 
We regard J as a polynomial in 8, and denote its leading coefficient* 
by J’(ay,...,B, 3). Then regard J’ as a polynomial in f,,_, and denote 
its leading coefficient by J”(a,,...,8,-2) and soon. Finally, we obtain 
a polynomial Jj(a,,...,8,) such that, if Jo(a{,..., 8) 4 0, the poly- 
nomial J(a,..., BO, By,...,B,) is not identically » zero. 

Ifa, ~ 0, 


7 


b b3 a > e27i SO, T+B 1,0. On T+ PJP — > ot y E2TESMLiyeoeyFyIIP 
© Be n 


ca Zn 
Since F, depends only on a4,...,«,, the number of solutions of 
x, J) F, = 0 is O(p™) by Lemma 2, and hence the number of sets 
(0),..+, &», 8) for which a, Jy Fj, 40 is p+»™— O(p™). It follows that 


f in +m __ O(pm)} py 7 ~ e27i S[f(x1,..-,L_)I/D - 
| ty Zn 
ais | (oY e2tiS [Hor 24 Bryn e+Bip | 
et kt 
aS oFe +0 
> E27 SN 12 +B ise On +B IP 1 
-E-BE[EE 2 
ii Toe ca J#0 
DY DY 0. Y e2PSM 01 H+ Bives ne + Bll? wl 
“5b 


By Holder’s maqnnty. this does not exceed 


92k- ‘2. 2 >| > a 2 E27 Sf, + Bj yee4y Xp L+By)I/P Pe 





SET eA ¥ 
| va 2k 
+. Q2k-1 QrriS[f(x, 7+B),..05 OX, 2+f,))/p 
2 — ter 12+ Bay iy t+By) 
ay An Bi &£ Be Bn 
aS oF +O J=0 


= 2, + &,, say. 
By Lemmas | and 2, 
= AF FEF. Foam) 
na "sage 
— O(pmlr+1)+2k(n—241-1/k))) — O(pmekn—2k+n—1))_ (4) 


y Hélder’s inequality, 


S, < 92k— Lgmin—1X2k—1) > a > > ae > 
Oy On 


D3 E27 S/(011 F+Biyooy Oy L+Py I/D ¥ 
snide Br 


oo 





* That is, the coefficient of the highest power of f,. 
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Suppose that 

f(x, x+B,,...,4,¢+8,) = Hvk+...tKat+F, 
and that J is the Jacobian of F,,...,F,, (lL <1 < ... < Tan < &). 
Then the number of solutions of the system of equations 

F,, = F,{a,...,B,,) = @,, (¢ = 1.,..,, 90) 
in the 2n unknowns 4,,...,8,, for which J + 0 is O(1), by Theorem 1. 
It follows that 
Ey = pmn WA DOS... S| emiStanehstarcailp "). 


ak @'2 


Hence, by the inequality from which we deduced Lemma 1, 


zy —_ O( pmn—DR2k—1)p2km) — O(pmekn—n +), (5) 
By (4) and (5) 
> bt 2 e27iS[f(ar = zip — O(pmekn—n +1)-(n +1)|/2k) ee O(pmra-ih)), 


Zr In 
the result required. 
I wish to express my hearty thanks to the referees for their kind 
suggestions. 














CONGRUENCE PROPERTIES OF RAMANUJAN’S 
FUNCTION 7z(n) 


By R. P. BAMBAH (Delhi), 8. CHOWLA (Lahore), 
H. GUPTA (Hoshiarpur), and D. B. LAHIRI (Calcutta) 


[Received 28 August 1946] 


1. FoLtLow1ne Ramanujan, we define 7(n) by the relation 
s r(n)x" = 2 II (1—a™)*4, 
where |z| < 1. : : 
It is known that (1) 
t(n) =0,,(n) (mod 691); (a) 
t(n) = no(n) (mod5), (b) 


where o;(n) is the sum of the kth powers of the divisors of n and 


a(n) = o,(n); 


zt(n) = 0 (mod 23) (c) 
if n=k (mod 23) and (5) = —1; 

t(n) = 0 (mod7) (d) 
if n = 0, 3,5,6 (mod7); 

o(n) = 1 (mod 2) (e) 


if and only if n is an odd square; 
t(jn—1) = o(jn—1) (mod)) (f) 
if j is a divisor of 24. 
Recently Lahiri has shown (2) how to determine the residues of r(n) 
modulo 3. His results are equivalent to the simple formulae 


t(n) =o(n) (mod3) if (n,3) = 1; (A) 
7(3n) = 0 (mod 3). (B) 


After reading an advance copy of Lahiri’s paper, we were able to 
cet new and simple proofs of the striking result (A), which, curiously 
enough, seems to have been overlooked by previous authors. 
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2. Proof of (A) 


Asin(3),wehave 1728 } 7(n)a” = G— 


where Q = 14+ 240 2 o,(n)a” 


R= 1—504 ¥ 0, a. 


and 


Now, if 4(2) be any power series with integral coefficients, then 


= « 


{h(a )jp— d( xP) = p 3 e, a”. 


T 
where eé, =90 (modp), p being a prime. 
In particular, when p = 3, we see that 
3 x 
P: o3(n)a”) = > ¢.2* 
1 
where c, =90 (mod3) if (n,3)=1; 
o3(U)o,(v)o,(w) = 0 (mod3) if (n,3) = 


OF = 5d, ". 


From (2), 
0 


1.€. 
Uurvct ,] 
u,v,wa1 


where, for n > 1, 
+2408 0,(U)o3(v)o,(w 


d 7200,(n)+3.240? S} o,(u)o,(v)+2 ) 
u+v=n u+v+w=n ~ 
U,v,wr1 (5) 


vt 
u,vrzl 


Again, by (3), we have 

p 03(U)og(v) = j29{07(n)—o,(n)}, 
Urv=n 

u,v. 


’(n,3) = 1, we have from (4), (5), (6) 
= 7200,(n)+ 1440{0,(n)—o(n)} 
2o,(n)} (mod 81). 


R? = 31,2" 
0 


so that, if 
(mod 81) 


d 


vt 


= 9{o,(n)— 
Moreover, from (3), 


1, 

1, = —10080,(n)+ 504? o;(u)o;(v) 
Ut+v=n 

u,v>l1 


where, for 7 


= —36o,(n) (mod81). 
Hence, from (1), (7), (8), if (n,3) = 1, then 


17287(n) = 9{o,(n)—20,(n)+-40,;(n)} (mod 81), 
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e., if (n,3) = 
3r(n) = o,(n)—20,(n)+40,(n) (mod 9), 

i.e. 3{7(n)—o(n)} = x {d?—2d’+4d°—3d} (mod 9). 
Now 
d— 2d? +-4d5'—3d = d(d?—1)(3+2d?—2d*) = 3d(d?—1)—2d*(d?—1)? 

= 0 (mod9). (10) 
Hence, from (9) and (10), if (n,3) = 

t(n) =o(n) (mod3). 


3. Proceeding on the same lines, we now prove 
7(n) =o(n) (mod8) 
a (0,2) <= ¢, 
Since d>=d*> (mod 8), 
we have Y a;(u)o;(v) = > a 3(u)os(v) (mod 8) 


u+v=n u+v=n 
u,v>1 u,w>1 


n)—o,(n)} (mod 8). 
From (5) and (6) for n > 1, we — 
d,, = 7200,(n)+ 1440{0,(n)—o,(n)} (mod 2°) 
= 208{20,(n)—o,(n)} (mod 2°). 


From (8) and (11) for n > 1, 


504? ‘ 
130 {o,(n)—o,()} (mod 2°) 


lL, = —1008o, ta) 5 ee 


= 1605(n) +55 {or(n)—o9(n)} (mod 2°). (13) 
5 
Now from (1), (12), (13), for n > 1, we have 
15.17287r(n) = 15. 208{20,(n)—o,(n)}— 15. 160,(n) —8{o,(n)—o4(n)} 
(mod 2°), 
i.e. 15.2167(n) = 11lo,(n)—300,(n)—5o,(n) (mod 2°), 
i.e. 24{a(n)—7(n)} = > {1ld?—30d'—5d*+-24d} (mod 2°). 


din 
Now, if (n,2) = 1, every divisor d of n is odd and 
11d?—30d>—5d3+4 24d = d(11d2—8)(d?—1)?—32(d?—1 
=0 (mod 2°). 
Hence, if (n,2) = 1, 
7(n) =o(n) (mod8). 
L 
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4. From §§ 2, 3 we have 


If j be a divisor of 24 and (n,j) = 1, then 
t(n) = o(n) (mody). 


It may be noted that (e) and (f) are particular cases of (D). 
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ON THE DETERMINATION OF THE JUMP OF A 
FUNCTION BY ITS FOURIER COEFFICIENTS 


By M. L. MISRA (Udaipur) 
[Received 19 October 1946] 


1. Let f(x) be integrable in the sense of Lebesgue over the interval 
(—7,7) and be periodic outside with a period 27. Let the Fourier 
series of f(x) be 


@ 
jay+ > (a, cosna+b, sin nz). (1.1) 
n=1 


Then the conjugate series of (1.1) is 
> (b,, cosnx—a, sin nx). (1.2) 


n=1 
Suppose that 


H(t) = Hflw+t)tfle—t)—2s}, 8 = a(z), 
We) =fee+)—fe—)—D, D= De), 
Lf og@\ YH au (ws 
f(t) = Tia) | (los) — du («> 0), 
i 
fo(t) = b(t). 


It is known* that 


I'(p 


Szasz}t gave the following theorem for the determination of the 
jump f(w+0)—f(#—0) or generalized jump D(x) of a function f(x) 
by its Fourier coefficients: 


t 


’ 7 6-4-3 
oa(t) = = = [ (lox) ff) du (B> a). (1.3) 
i 


THEOREM 1. Jf there exists a number D = D(x) such that 
t 


| ¥(t) dt = o(t) and [ w(t)| dt = 0, 
0 


0 


- _ 9 
as t > 0, then lim {8}, (x) —S3(x)} = —=— D(z), 
n> 7 
where S1(a) is the sequence of the arithmetic means of the partial sums 
of the conjugate series. 


* Wang (6). t+ Szasz (5). 
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Theorem 1 was generalized by Chow* in the form of 
THEOREM 2. Under the same hypothesis as in Szdsz’s theorem 


lim (53, (x) —Sa(x)} = “°F? Dixy, 

no 7 
for « > 0, where S%(x) is the n-th Cesdro mean of order « of the conjugate 
series. 

Further generalizations of this theorem have been given by 
Minakshisundaram.} 

Definition. A series > c, is said to be summable by Riesz’s 
logarithmic mean of order k > 0, or summable (R, logn, k) to the 
sum s, provided that 

i. a oo > (lose 
e() (log w)* 4 "“— 
tends to a limit s as w > 00. 
The object of this paper is to prove the theorems: 


THEOREM A. If 


as t > 0, then lim [R,(2w)—R,(w)] = me Su D(z), 


wo Qa 
where R,(w) is Riesz’s logarithmic mean of order unity of the conjugate 
serves. 


THEOREM B. If 


j pla allen) 
t 


as t > 0, then 
— Ps D ! ; 
lim [Fig 1(Aw)— Fg a()] = G42) 8" (A>1; B>0), 
where Rg,,(w) is Riesz’s logarithmic mean of order (B+-1) of the con- 
jugate series. 

Theorem A is the case A = 2, 8 = 0 of Theorem B. 

I also prove the following generalization of the well-known 


* Chow (1). + Minakshisundaram (4). 
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theorem* on the order of magnitude of the partial sums s,,(x), &,(x) 
of the Fourier series of f(x) and its conjugate series respectively: 


THEOREM C. (a) If 


as t > 0, then 


(b) If 


as t > 0, then 


[t is easily seen that, if 


as t > 0, then 


On the other hand, if 


fla = o( lox) 
t vl i 


t 
t 


1 
then [ \p(t)| dt = o(tios 
0 
Hence Theorem C generalizes the theorem on the orders of magni- 
tude referred to above which holds under the conditions 


t 
| Isl de = of), 


0 


2. I shall make use of functions S,(t) given by 


1 
S,(t) = | (vos; J sintudu (k > —1). 
0 


* Hardy and Rogosinski (3) (49, Theorem 64). 
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It is known* that 


© [08,0 )] = kS,1() (k > 9), (2.1) 


Y P(r+s+2) [s ] P 
A = - S.(ut ] _ 1 Pr. —1). 
, (2.2) 
We shall require the lemma 
Lemma. If k > 1, then forall w > 1 
k 


w log w 


Ri(w) = [R,,-1(w)—B;,(w)]. 


We have 


s d log n\* 
Fo) = J] >. (ing =) | 
k log n\*-1 
oe Pe *) (log n)ec,, 


log WwW 
' oa \e-2 k 
- k >((- log n - ae log n c., 
w log w ~,\ log w log w 
which proves the lemma. 


3. Proof of Theorem C 
(a) We have 


exe =i n — 5 
8, mo s = L 
sin ie 

7 / In 7 


aT Ih (a ae +0(1) 


(+ I,)+-0(1), say. 


7m 


| <n { |p| ae 


0 


= | fe al + f| { \g 


i 0 0 ‘t 
=o(logn) (no), 


* Wang (6). 
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(6)| dt = o(loe:) ast>0. Also 


i,= of | 40a) = 0(logn). 


Hence 8, (x) = s+o(logn)+0(1) 
=o(logn) (no), 
which proves the part (a) of Theorem C. 


(6) We have 
= | { f(x+t)—f(a—t)} cot }4(1—cos nt) dt +-0 (1) 
0 


7 
) 
= [ cot 44(1—cos nt) dt + 
2a = 
0 ; 
Tin 7 


+n | + [ }yeeot4(1—cos nt) dt +0(1) 
: 4 


m/n 


J,+J,+J3+0(1), say. 
Now* J, = Pog 2n-+-y]+0 (1), 
7 


where y is Euler’s constant. Also 
mn 
Jy) <n | |p(t)| de 
0 
= o(logn) (no), 


if [ ty dt = o(lo:), as t— 0. 


i 


Again J, = o( | Ha) = 0(logn). 


Hence 


&,(x) = 7 tog n-+-0(1)-+0 (log n) 


D 
~ —logn, 
7 


which proves the part (b) of the Theorem C. 
* Hardy and Rogosinski (3), 47. 
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4. Proof of Theorem A 
We have, if 


<w 


> (6, cosna—a, sinnxz) = > B,, 
nw 


7 
7 


= - | { f (x-+-t)—f(a—t)}cot (1 —cos wt) dt +0 (1) 
27 . 


0 


Df lf 

= cot 3¢(1—cos wt) dt +-— | ¥eeot 4t(1—cos wt) dt +0(1) 
a7 2a 

0 0 


D 1 f 2 
— (log 2w+y)+5 [ 40 (eo jt— 7) 1—cos at) dt dt + 
7 oT 


Y 
0 


t 
+2 [oy 2" at +0(1) 
7 t 
0 


) 
1—cos wt 


dt +0(1) 


P og w +4 [—(1—cos wt, (0) +2 | i, (t)sin wt dt +-O(1) 
7 7 7 
0 


a OE | is, (t)sin wt dt +O(1) (w+); 
7 7 
0 


wW 


] = dx 
| Raj 
paw 


= [ °8 ’ dx oe [ ax fv t)sin xt dt +O(1) 
mlogw) 2 stag | 
1 


7 @w 


B 
— log w + 


a ~*~ + apes J #0 t) dt fanmse+00 


= Jog w+ f a0 1— oe at +0(1), 


7 ice w 
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7 
1—cost 
| 28 at = oc, 
0 


,(t) being integrable (L) in (0, z). 
Hence, by the lemma, 
= a - 

7 Ry(w) = —7__ [R,(w)—R,(w)] 


dw wW log wW 


—: 
~ wlogw|2z 


i} te +0] 
or 
D ] 1 ] 
= --+——_ [~+ ——__J —}. 
2w  wlogw w(log w)? a 
Now I = o(logw), 
as in the proof of Theorem C. Also 


w 


JI < | fo yb, (t) ae dt +| +] [sn Sete 


my 
lelao 


a7/w 


=w — 0(log 7) dt +- [ 7 (toe) dt 
“t Jf t 
0 7/w 
0 (log w)-+-of (log w)?} = 0 {(log w)? 3 
Hence 


Now 


which proves Theorem A. 
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5. Proof of Theorem B 
We have 


—logw +" | i), (t) sin wt dt +0(1). 
“4 


Also, for k 


w 


w\* k [[,_ o)\**R, (2x) 
fe) = a) Bos a let a 
x() (log w)* > (ie °] B, (log w)* | ( ™ 4 a. 
i 


now 


DD: & r -Llog x 
<i zm (log w)* . (los =)" = - dx + 
1 
k r k-1 . , 
m(log w)* ‘ (los) da | ,(t)sin xt dt + O(1) 
1 0 


oa Diogw | k ub, (t) dt | (es: “\"" sin xt dx + O(1) 


(k-+-1)a Alo oF | 


Dlog w k eee : 
a (k+ 17 | a(log a) :) sin at dx + O(1) 


D log w 
_ Piogw , % t)S.._,(wt) dt +-O(1). 
(EFI) t aflogaye J Pi Oealed) a +00) 
Hence, by the lemma, we have 


Bo ,s(o) = PE) (Bw) Rp (co) 


w log w 


D 


co 
w log w 
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Now, integrating by parts h times where h = [f], if 8 is not an 
integer, we have 


7(log w)P+1 
—BE+1) 


“a | why (t)Sp_y (ot) at 
0 


= [—tp,(t)Sp_s(wt) |} +(B—1) j o(t)Sp_o(wt) dt 
since tyb,(t) > 0 as t > 0*] 


= (B—1) | do(t)Spa(wt) dt = 
0 


7 


= (8—1)(B—2)...(B—h) J Yrsalt)Sp_p—a(wt) dt 


lid r id 
_ (B—1)(B—2)...(B—h) i Sp_p-a(wt) dt | (los) dale) 
3 é 


T'(h+1—B) 


3G -B 
_ ee” | vt (oe ‘) Sp_p-(wt) dt 


=F) fo u)S,(cou) oe COS we 2 


by (1.3) and (2.2). If 8 is a positive integer, we obtain the same 
formula after integrating by parts B—1 times. Hence, if (1.5) is 
satisfied, we have, as in the proof of Theorem A, 


Ru o(*) (w > 0). 


Similarly, by (1.5), 


ee ag a= fa a(t )So( (cut) 


- = o( j \dpaa(t)|-—<° i) = o(2). 


* See Hardy (2). 
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pu D 1 
Therefore Rpis(w) = @42)ne +o (=); 


and 


Aw 
Rg.,(Aw)—Rg,1(w) = | Rgsr(t) dt 


D 
= Grae erro », 


which proves the Theorem B. 


I am much indebted to Dr. B. N. Prasad for his kind interest and 
advice in the preparation of this paper. 
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WELL-POISED BASIC HYPERGEOMETRIC SERIES 


By W. N. BAILEY (London) 


[Received 11 September 1946] 
1. Introduction 
DovuGALL’s theorem for the sum of a well-poised hypergeometric 
series ,F, of the ordinary type was given in 1907. In 1932, I gave 
in this Journal* the formula 
pl® 1+44, b, ¢, d, é, f; 
, | ja, 1+a—b, 1+a—c, 1+a—d, 1+a—e, 1+a— 
a l(1+a—c)(1+a—d)'(1+a—e)'(1+a—f) - 
~ Td+a)P(o—a)P(1+a—d—e)P'(1+a—c—e)I'(1+a—c—d) 
. P'(b+c—a)l'(b+d—a)P(b+e—a)l(b+f—a) _ 
l(i+a—c—f)C(1+a—d—f)l(1+a—e—f) 
_ TU+2b—a)P'(b6+c—a)P(6+d—a)P'(b+e—a)P(b6+f—a) 
N(1+6—c)P(1+6—d)(1+6—e)(1+6—f) 
l'(a—b)T(1+a—c)'(1+a—d)l(1+a—e)l(1+a—f) = 
P(b—a)P(1+a)P(e)P(a)P(e)l(f) 
E 2b—a, 1+b—}a, b, b+c—a, b+d—a, 
on b—ja, 1+b—a, 1+b—c, 1+6—d, 
b+e—a, b+f—a; a: 
1+b—e, 1+6-—f 














where 1+ 2a = b+c+d-+e-f. 
When one of the parameters c, d, e, f is a negative integer, this 
reduces to Dougall’s theorem. 
In the same paper I obtained a formula} connecting four well- 
poised series of the type ,F,. If we write 
P(1+a)P(6)P(e)PA@)Te)l(AL@)T(A) 
hl+a—b)r(l+a—c)[(1+a—d)['(1+a—e) x 
x P(1+-a—f)P(1+-a—g)P(1+-a—h) 
1+ 3a, b, C, d, e, 
ja, 1+a—b, 1+a—c, 1+a—d, 1+a—e, 
f 9; h; 
l+a—f, 1+a—g, 1+a—h]’ 


* Bailey, 2 (7.4). See also Bailey, 1, § 6.6 (3). 
+ Bailey, 2, § 9. For the form given here see Bailey, 1, § 6.8 (3). 





V(a;b,c,d,e,f,g,h) = 


x ml” 
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the formula is 





cosec(b—a)n{V (a; b,c, d,e,f,g,h)— 
—V(2b—a; b, b+c—a, b+-d—a, b+-e—a, b+f—a, b+g—a, b+h—a)} 


P(e) (d)l(e)l'(f+b—a)l(g+b—a) 
x Pea ay ae 
P(h-+b—a)cosec(1+-a—f —g—A 


P(1+a—g—A)P'(1+a—f—h)P(1+a—f—g) 
x {V(1+2a—c—d—e; b, 1+a—d—e, 14+a—c—e, 14+a—c—d, f, g, h)— 





—V(2b—2a—1+c+d-+e;b,b—a-+c, b—a+d,b—a-+e, 
1+a—g—h, 1+a—f—h,1+a—f—g)}, (1.2) 
provided that 2+3a = b+c+d+e+ft+g+a. (1.3) 
By iterating this formula, I showed* that 
cosec(b—a)n{V (a; b,c, d,e, f,g,h)— 
—V(2b—a; b,b+c—a,b+d—a, b+e—a,b+f—a,b+g—a,b+h—a)} 
+ P(e) P(d) P(e) P'(f )P(c+b—a)P(d+b—a) 
~ T(l+a—c—g)P'(1+-a—d—g)'(1+a—e—g)I'(1+a—f— 9) 
T'(e+b—a)(f+b—a)cosec(g—h)x 
5 I'(1+a—c —A)fi+a-—-d— h)V(1+a—e—h)T( (1--a—f—h) * 
x {V (b—g+h; b, 1+-a—c—g, 1+-a—d—g, 1+a—e—g, 
1+a—f—g,h,h+b—a)— 
—V(b+g—h; b, 1+a—c—h, 1+a—d—h, 1+a—e—h, 
1+a—f—h,g,g+b—a)}, (1.4) 
provided that the condition (1.3) is satisfied. 
In (1.2) and (1.4) any one of the series ,F, is of general well-poised 
type, except for the second parameter and the restriction that the 
sum of the denominator parameters exceeds that of the numerator 


parameters by two. When h is a negative integer, these formulae 
reduce to formulae connecting two terminating well-poised series 


oFy-t 


It is now twenty-five years since Jacksont{ gave the basic analogue 


* Bailey, 1, § 7.6 (2). + Bailey, 1, § 4.3 (7) and § 7.6 (1). 
t Jackson, 4. See also Bailey, 1, § 8.3 (1). 








| 
| 
| 
| 
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of Dougall’s theorem. In 1929 Watson gave a transformation of a 
terminating well-poised ,®, into a Saalschiitzian ,®;, corresponding 
to Whipple’s transformation of a well-poised ,F,, and in the same 
year I gave a transformation* of a terminating well-poised ,,9, 
which is exactly analogous to the particular case of (1.2) when h is 
a negative integer, and gives a generalizationt of Watson’s trans- 
formation of an ,®, vshen the series does not terminate. 

It has thus been obvious since 1932 that formulae corresponding 
to (1.1), (1.2), and (1.4) are probably true for basic series. All these 
formulae, however, were discovered by means of contour integrals 
of Barnes’s type, and it did not appear obvious how to attack the 
problem for basic series. I now find that the basic analogue of (1.1) 
is an immediate result of a formula which I published in this Journal 
ten years ago.{ 

In the present paper I give the basic analogues of (1.1), (1.2), and 
(1.4). The proof of the analogue of (1.2) was suggested by the method 
of obtaining transformations of integrals of Barnes’s type§ by which 
(1.2) was discovered. This proof, when written out for (1.2) itself, 
consists in using the series equivalent to the integral of Barnes’s type, 
but the details both for (1.2) itself and its analogue were found to 
be surprisingly troublesome. 


2. Notation 
In this paper I use the notation 
(a),, = (1—a)(1—aq)(1—aq?)...(1—ag"“"), (@)o = 1, 
where |q| < 1. 
[ also use the function II(a) defined by 


1/IIl(a) = 


n 


(1—aq"). 


8 
it 


3. The basic analogue of (1.1) 

If 
: res Il (aq) ig st 
xl4;6,¢,d,¢.f] = rao7pyit(aq/e)ll(ag/d)ll(aq/e)l1(aq fy 1 (a®g bodef) * 


x .o,|% qva, —qva, 6b, ec, .-  €& 720 a°g?/bedef 
alii va, —va, ag/b, ag/c, ag/d, ag/e, aq/f 4 





* See Bailey, 1, § 8.5 (1). + Bailey, 1, § 8.5 (3). 
t Bailey, 3. § See Bailey, 1, § 6.7, or 2, § 8. 
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then it is known that* 








xl@;6,c,d,e,f] 
TI (ag/def )1 (def/a)T1 (6d /a)T1 (be/a)T1(g/e) M1 (bf/a) 
Le _xlefle;e, f, aq/be, aq/ed, ef/a] 5 





~~ II(aq/b) 1 (6/a) 1 (ag/ef 0 (aq/df )I1( (aq/de) Il (a2g2/bedef ) 

‘: bx[b?/a; b, be/a, bd/a, be/a, a EERE 

all (aq/bc) 11 (d)I1 (e) 11 (f M1 (a*g/bdef ) 1 (bdef/a*)’ 

If we take a*q = bedef, the first series ,P,-on the right-hand side 
of (3.1) reduces to a ,®; which can be summed by the formulat 


o,| ™ b, iy d; agibea 
va, —va, ag/b, aq/c, aq/d, 


(3.1) 


__ II (a@q/b) I (aq/c) 1 (aq/d) 1 (ag/bed) (3.2) 
~~ Ti (aq) 11 (ag/be) 1 (ag/bd) 11 (ag/ed) * ii 





We thus find that 
o.|% qva, —qva, 6b, C, d, ee 
sie va, —va, ag/b, ag/c, aq/d, aq/e, aq/f q 
II (ag/c) 11 (ag/d) 1 (ag/e) 11 (ag/f) 11 (be/a) 1 (bd /a) 11 (be/a) 11 (6f/a) 
™ Taq) (b/a)I1(aq/de)M(aq/ce)I1(aq/cd) 1 (aq/ef )I(aq/df M1 (aq/ef) — 
I] (67g /a) 11 (be/a) 11 (bd/a) 11 (be/a) 1 (bf/a) x 
x I (aq/c) 11 (ag/d) 1 (ag/e) 1 (ag/f )11 a 








~~ TI(bq/c) 11 (6q/d) 11 (bq/e) 1 (bq/f)T1(aq) 11 (c) 11 (d) 1 (e) 11 ( f)T1(6/a) 
a ‘ee qb/Na, —qb/va, b,  be/ja, bd/a, be/a, bf/a; 
nla | b/va, —b/va, bq/a, bq/c, bg/d, bg/e, bq/f 


where a2q = bedef. (3.3) 


This is the analogue of (1.1). 


4. An alternative form of (3.3) 
Since II(aq”) = (a),,II(a), we can write (3.3) in the form 


IT (ag™) I (g"*4va)I1(— q”*1va) 


(Dll (q"va)I1(— ga) II (ag, b) 


ad aes IT (cg™) II (dg) II (eq™) 11 ( fg™)q” 
TI (ag™**/c) 11 (ag™*1/d) 11 (ag™*1/e)T (ag”*1/f ) 


m=0 





* Bailey, 3 (5.1). 
+ This is the particular case of Bailey, 1, § 8.5 (3), when ag = de. 
{ Corresponding to Bailey, 1, § 6.6 (2). 
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_ To) re )I(e) Hen IT e/a) El (bd fa) 1 (be/a) Ioffe) 
II (ag/de) 11 (ag/ce) 1 (ag/ed) 11 (ag/ef) x 
x I (ag/df) Il (aq/ef )11(b/a)11(aq/b) 
4b + Nth ee te Sel 
ntl (bq”/ Va) 11 (—bg™/Va) Tl (bg *4/a) 


y, M(begn/a) II aaran (beg (a) NOfa"/aya™ (4.1) 
II (bg **/e) TI (bg **/d) 11 (bg *4/e) TI (bg *4/f)” 














Now replace a, c, d, e, k, b by k, ke/a, kd/a, ke/a, a, at, where, as 
usual, k = a*q/cde, and we get* 
Ii (at) 11 (ct) 1 (dt) 11 (et) 11(1/t) 
Il (atg/c) 11 (atg/d) 11 (atg/e) 11 (at/k) 1 (kg/at) 
II (c) II (d) 11 (e) 


~~ [(ke/a)I1 (kd /a) 11 (ke /a) 11 (a/k) 
x[> of q™ Il (kq™) I(g"*1vk) Tl (—q m+1/i)TT (atq™) 














mil (g@vk) Tl (—q@ vk) (kg™*1/at) Tl (ag™*1/c) 


I(keq”/a)T1 (kdg/a) 11 (keg”/a)T1(q"/t) 
(ag *¥/d) (ag /e) 11 (tg) 


_at = ql ( at qk) saat —atg”*1/vk) 
ED, 


=, ' 








Il (atg/Vk)I ee at? gm +) 
Il (ag /k) TI (ctg™) II (dtg) 11( etq™) II (atg™ i (4.2) 


Il (detq™) 11 (cetg™) Tl (edtq™) 1 (atg™*1/k) 
If t = q", where n is a Lape integer, this reduces to 
II (ag”)I1(cq") 11 (dq”) To, 





I (ag"* +1, Siena EM aq" /k) 11 (kq'- n/a) 
II (c)II(d) 11 (e) 


"" _—_ (kd/a)TI(ke/a)T1(a/k) 


q™ ii (kg™) 1 (q¢ m+1/k)T1(—q™*1vk) II (ag" +") 
(Vm El (gq vk) TT (—g vk) TT (kg—" + /a) I (ag™*4/c) 














eee (kdgn/a) IT (keg”/a)(q")m (4.3) 
(ag™*1/d)T1( aq™*1/e)IT (kqgm™+"+1) é ' 


This formula is, of course, pel to Jackson’s analogue of 
Dougall’s theorem. 


j 


* Corresponding to Bailey, 1, § 6.7. 
3695-18 M 
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5. A transformation from (4.3) 
Using the identities 
T(Ag-) = (—1)g'™*11(A)/A?(q/A),, 
I1(4q")I1(q'-"/A) = T1(A)I1(q/A)(—1)"Anginen®, 
we get from (4.3) 
an II (agq”) 11 (eq") IT (dg”) 11 (e gM (oy 9") ee q” a” 


I (ag"**/c) 11 (ag"*4/d) TI (ag"**/e) IT (o, g)T(aq"*4/b) 
Te)Ma)I I(e)I1(kq/a) 


= Ti (ke/a) Il (ed /a)L1(ke/a) 
n ql sn m+1vk)I1(—q q” vk) II (ag"*™) Il (keg™/a) 
(Vil (g@vk) 1 (—q" vk) Tl (kg”"—"*4/a) TI (ag™**/c) 
I (kdg™/a) 11 (keq/ a)(q~ es 
* Ti(ag*1/d) 11 ( aq™* L/e) I (kg™t"+1) 
IT(p, 9") (p2q ,q" a (—1y(7)" aime 
(q)n l(a, 9") 1 (aq"*1/b) k 


Writing n = m+», and using the identity 


— 1)” "(Q) m+p 





n=0 n=0 


—m—p\ — 
(q Ss a (q \»p qim@p+m+1)? 


the right-hand side becomes 
He) (@)i(e) 
aia aon i 
— (kg) 1g m+1 vk)II (—q"*1vk)I I (aqg?”) 11 (keg /a) 1 (kdg™/a) 


(2) mn Th(g@vk)l (—q"vk)I T(ag™*t/e II(ag”* +1 /d) 


) 
I (keg™/a) IT (p, q")Il (p. q”) __ ax m ” 
I (ag *¥/e) 11 (kg?"*4) TI (o, g™) I (ag™*4/b) \& 
aq?", a [k, Py g”, Po q” : 
@ : ', ‘ 
kq?m+1, o} ¢, aqg™*1/b | 
In this result there may be any number of the symbols p, o. In 
particular we find that 
oe II( (aq") )IT( q"*1Wa) Il (—q"*+1 va) IT (bg") TI (eq”) I (dq”) 
>i (7), 11 (q" va) 11 (—g" va) I (ag"*1/b) I (ag"*1/c) 1 (ag"*1/d) 
ss II (eq”)IL( fq") U1 (gq" IL (Ag”) 
I (aqg”*+*/e)I1(ag"**/f II (aq” +1/g) IT (aq"*1/h) 





q” 
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Ei II(c)I1(d) 11 (e) 
Il (ke/a) 11 (kd/a) 11 (ke/a) * 





TI (kg) 11 (g++ vk) Tl (—q™*1vk) x 
S x TI (ag? +1) TT (keg™/a) I (kdq™/a) 
Dm (g" vk) I (—gq" vk) I (ag™**/c) 1 (ag™**/d) 





m=0 


- II (keg™/a)I1(bg™)I1( fq™) II (gq™) 11 (hg™)(aq/k)™ 
II (ag™**/e) I] (ag™**/b) 1 (ag *1/f 11 (ag */g) Il (ag™**/h) 





© aq”, qg™* va, —q™*1va, bq”, fa”. 
X 7 mala —g™Va m+1/h m+1/f 
qd va, q ? aq ais | aq ? 


99" hg”, ajk; 
aq” *1/q, ag” **/h, kg? (5.1) 


The series on the left of (5.1) is, of course, a multiple of a well- 
poised series ,,P,, with numerator parameters a, gva, —qva, b, c, d, 
e, f, g, h. The change in the order of summation is easily justified 
because of the powers of q that occur. 


6. A transformation from (4.2) 


Now in (4.2) take ¢ = bq"/a, and we obtain the sum 








=. 11 (bg")I1(beq” /a) TI (bdg”/a) 1 (beg”/a)I1(p, bg"/a)T1(p,bq"/a), , 
Da, TI (6qg"*1/c) 1 (6g */d) 11 (bg”*4/e) TI (bg *4/a) I (a, bg”/a) 


as the sum of two double series. In the first we simply change the 
order of summation and obtain 


IL (c) 11 (d) 11 (e) 11 (6/k) 11 (gk/b) 


Il (ke/a) 11 (kd /a) 11 (ke/a) 11 (a/k) 11 (bg/a) 11 (a/b) 





aad ql (kg) Il fymtty k)T ( —gmtly kj) (bq™) TI (keq™/a) 
(9) ll (g’vk)Il(—q@vk) IT (kg +, b)Il(ag™**/c) 





m=0 


Il (kdg”™ /a) 11 (keg™/a) II (ag™/b) IT (p, 6/a) 11 (p, b/a) 
Il (ag”*1/d) I (ag *4/e) Il (kbg™*1/a)I1 (oc, 6/a) 


oe" Sa-“/k, mdi, pablo; , 
bq'-™/a, kbg”*1/a, o,b/a p 





x 
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In the second double series on the right we put m = p—n and then 
change the order of summation, and we obtain 


b (c) IT (d) I (e) 1 (6/k) 11 (gk/b) 
—k eaves a)I1 (ke/a)I1(a/k)I1(bg/a) I (a, ib) * 


«> q? I1(b?g”/k) IT (bq?**) vk) IL (—bgq? +4/Vk) 11 (ag? /k) 


(7), 11 (bq?/Vk) TN (—bq” vk) 11 (b?g”*/a) 
Il (bcq” a) IT (bdg? /a) 11 (beq” a)I1(bg”) Il (p, b/a)I1(p, b/a) 


Il (bdeq” /a) Tl (beeg” a) I (bedg” /a) 11 (bg? +4/a) II (, b/a) 
p,b/a, pob/a, b*qP/k, qr; 
xo] o,b/a, kq'-¥/a, b?g?*/a 

Again there may be any number of the symbols p, o. In particular 
we have 

II (6°q” a)Il(q"**b, va) II (— gq’ 1b va) IT (bg” IT (beq”, a) 
> (q), 1 (q"b/Va) I (—q"b/Va)T (bg”*1/a)¥" bg" *1/c) 

Il (bdq” /a) 11 (beq” /a) I (bfq"/a) 11 (bgg”, a) 11 (bhg"/a)q” 

Il (bg”+1/d) 11 (bg”*1/e) IT (bg”+4/f ) I (bg”+4/g) Il (bg"*+4/h) 
__ I (c)I1(d) 11 (e) 1 (b/k) I (gk/b) 11 (6°g/a) I (bf/a) 11 (bg/a) 11 (bh/a) 


~ W(ke/a) I (kd/a) 1 (ke/a) 1 (a/k) 1 (a/b) 11 (bg /f )T1( at (bq/h) 


gq" il( kq) IT(g m m+ivk)T(—g"*1vk) IT (bg) bey 
«(2 (q),, 11 (q"vk) Tl (—q" vk) I (kg *2/b) 1 (ag"*1/c) 
Il (keg™/a) 11 (kdg™, a) II (keg™/a)II(ag™ /6) 
Th (ag™*t/d) I (ag@*1/e) Tl (kbg"*4/a) 


n=0 


m 


wails ‘b?/a, gb/va, —qb/va, bq”, bq-™/k, 
es bva, —bva, bq'-™/a, kbq™/a, 
bf/a, bg/a, bh/a; 
bg/f, bg/g, bg/h 
q? 11 (62g? /k) TL (bg? **/vk) 11 (—bq?**/ Vk) IL (ag? /k) 
a >) (q),, 11 (bq? /Wk) IL (— bq” / vie) 11 (b?g?+4/a) ; 
II (beg? /a) 11 (bdg? /a) 11 (beg? /a) 11 (bq”) 
Il (bdeq” /a) 11 (beeg”/a) 11 (bedq” /a) Tl (bg? */a) 
— b?/a, qb/Na, —gqb/va, bf/a, bg/a, bh/a, 
os b/Na, —b/va, bq/f, b/g, bg/h, 


bigrlk,  q7?; 
kq-/a, b2g?ia al}. (6.1) 
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7. The basic analogue of (1.2) 
If a*q* = bedefgh, (7.1) 
the second series ,®, on the right of (6.1) can be summed by Jackson’s 
analogue of Dougall’s theorem. Also, under this condition, the first 


series «D, on the right of (6.1) can be combined with the series ,, on 
the right of (5.1) so that we can use (3.3) with a, b, c, d, e, f replaced 


by aq”, bq”, fa”. gq", hq”, a/k. 
Writing 
W(a;b,c,d,e,f,g,h 
II (aq aii (c) IT (d) IT (e)I1(f )IL(g) IL (A) 
TI(ag/b)I1 (ag/e) 11 (ag/d) 1 (age) (aq/f )I (aq/g)T1(agjh) * 
xX sao a, qva, —qvs, 6, 6& €@ « Jf, gg, BF 
iB. va, —va, ag/b, aq/c, aq/d, ag/e, aq/f, aq/g, ag/h, 4 
we find, after considerable reduction, that 
I1(b/a)I1(ga/b)[ W (a; b, c,d, e, f,g,h)— 
— (b/a) W (b?/a; b, be/a, bd/a, be/a, bf/a, bg/a, bh/a)) 
II (c) IT (d) 11 (e) 1 (bf/a) 11 (6g/a) 11 (bh/a) 11 (6/k) 11 (gk/b) 


ie ~~ TI(ke/a) {1 (kd/a) 11 (ke/a) 1 (aq/gh) Il (aq/fh) 11 (aq/fg) 
x [W(k; b, ke/a, kd/a, ke/a, f, 9, h) 
— (b/k) W (b?/k; b, be/a, bd/a, be/a, bf/k, bg/k, bh/k)|, (7.2) 








where k = a*q/cde, and a*q? = bedefgh. 
This is the analogue of (1.2). 


8. The basic analogue of (1.4) 
By iterating (7.2) we obtain a formula equivalent to 
I1(b/a)I1(ga/b)[ W (a; b, c,d, e, f,g,h)— 
—(b/a) W (b?/a; b, be/a, bd/a, be/a, bf/a, bg/a, bh/a)] 
II (c) IT (d) IL (e) I1(f) ec /a)I1(6f/a) 
I] 


II (ag/eg)T1 (aq/dg)T1 (ag/eg)T1 (agq/fg) 11 (aq/ch) x 
x Il (aq/dh) Il (ag/eh) I (aq/fh) 


x I (g/h) I (gh/g)[ W (bh/g; 6, ag/eg, aq/dg, aq/eg, aq/fg, bh/a, h)— 
—(g/h) W (bg/h; 6, aq/ch, aq/dh, aq/eh, aq/fh,bg/a,g)], (8.1) 


where again a*g? = bedefgh. This is the analogue of (1.4). 
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9. Conclusion 

From (7.2) and (8.1), when one of the parameters is of the form 
q~", where v is a positive integer, we can deduce relations connecting 
two or three well-poised series of the type ,M,. By letting n > co 
we obtain relations connecting two or three series of the type ,®, 
in which the parameters are unrestricted, apart from convergence 
conditions and the special forms of the second and third numerator 
parameters. My transformations of terminating ,,®, mentioned in 
§1 are particular cases of (7.2) and (8.1). One of these contains 
Jackson’s analogue of Dougall’s theorem as an obvious particular 
case, just as Dougall’s theorem itself can be deduced from (1.2). 
Thus (7.2) and (8.1) include as particular or limiting cases all the 
known results for well-poised basic series. 
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THE TWO-DIMENSIONAL AEROFOIL IN A 
BOUNDED STREAM 


By A. E. GREEN (Durham) 
[Received 17 December 1946] 


1. THE two-dimensional perfect-fluid theory for finding the lift and 
moment forces acting on an aerofoil in an unlimited fluid is well 
known for a wide variety of aerofoil shapes. Very few exact solutions 
exist when the fluid is limited by one or more boundaries or when 
the aerofoil is near other bodies. The forces acting on a flat-plate 
aerofoil which is in a stream bounded by one or two straight parallel 
walls have been found by Tomotika and others (5-12) using con- 
formal transformations. The forces on a circular-are aerofoil in the 
presence of a plane wall were found by conformal transformations 
by the present writer (1). Tomotika has shown that a series form 
of solution for the flat-plate problems may be deduced from his exact 
analysis, the series being very suitable for numerical computation 
provided that the plate is not too near the wall. The flat-plate 
problems have been re-examined by Havelock (2) by quite a different 
method which leads directly to expansions for the lift and moment 
forces acting on the plate. His analysis actually applies to an elliptical 
aerofoil of which the flat plate is a limiting case. The method used, 
however, appears to be somewhat complicated, and it is of interest 
to re-examine the problem once more since the series solution can 
be obtained fairly easily by using quite elementary ideas. Moreover, 
the present method applies not only to the elliptical and flat-plate 
aerofoils but also to the usual more general types of aerofoils. In § 2, 
therefore, the problem of the fairly general aerofoil in the presence 
of a plane wall is discussed, and the extension to an aerofoil in a 
tunnel bounded by two parallel walls is indicated in § 7. 

In discussing the effect of the ground on the lift of a monoplane 
aerofoil Tomotika (10) found that his results for a flat plate were in 
qualitative agreement with experiment. The results for a circular-arc 
aerofoil differed considerably in some cases from those for a flat-plate 
and from experiment (1). The work of § 2 when applied to the circular- 
are aerofoil shows more clearly why the flat-plate and circular-arc 
aerofoils are different, and it also enables us to see how the thickness 
of the aerofoil affects the results. Only three terms of the series 
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expansion for the lift acting on the general aerofoil in the presence 
of a plane wall are given in the present paper so that results obtained 
from these terms must be accepted with caution. It appears, how- 
ever, that, for large distances of the aerofoil from the ground, camber 
and thickness both add to the decrease in the lift which is found on 
a flat plate compared with the lift in an unbounded stream. As the 
flat plate gets much nearer the ground it is known that, for small 
angles of attack, the lift eventually increases above the value found 
in an unbounded stream. In this case the present work indicates 
that the effects of camber and thickness on the flat-plate results 
appear to act mainly in opposite directions so that we might expect 
that the flat plate would give a reasonable approximation to actual 
aerofoil shapes. 

Havelock (2) also obtained expansions for the forces on a flat plate 
when the stream is bounded by parallel free surfaces, or by one rigid 
and one free surface, by taking the boundary condition in an approxi- 
mate form. The methods of the present paper could be used for 
similar problems when the aerofoil has a more general shape. 

In §7 the method of solution is indicated for the problem of an 
infinite cascade of similar aerofoils in a uniform stream, a problem 
of some interest in airscrew theory. The exact solution is known for 
a cascade of flat-plate aerofoils but not for other shapes (13). 


2. Consider the two-dimensional motion due to a cylinder whose 
cross-section is the curve y = 0, € = 0 ~ —2z in the transformation 


—— e~tt 3 ay eint (¢ _— £+%n) 
n=0 


placed in a uniform stream bounded by one plane wall y = —b. 
The origin of coordinates is inside the cross-section C of the cylinder. 
If c is the stream velocity parallel to the negative direction of the 
x-axis and if a circulation K = 27A, also exists round the cylinder, 
a suitable complex potential function w which makes the boundary 
wall a stream-line is 


w= cz+iA log z—log(z -++-2tb)}+- > {A, + ose) 
gn (z 2a 


n=1 


where a bar placed over a quantity denotes the complex conjugate 
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of that quantity. Apart from a constant, equation (2.2) may be 
written in the form 


w = cz+iA,logz+ , (2.3) 


where = (2.4) 


grt 7 
and "1. = Eze - ‘ " (2.5) 
Using (2.3) in the formulae 


(2.6) 


M= —IpR b= (=) dz, (2.7) 


for the forces (X,Y) and the couple M acting on the cylinder, we 
find that 


Y+:1X = 2mp\cAy+Ag B+ iS m(n+1) A, Bru}; (2.8) 


M = 2npRi{eA,+ > nA, B,,| (2.9) 
n=l 


nj° 


3. It is now necessary to satisfy the boundary condition on the 
aerofoil, Using (2.1) we may express 2~’, log z, 2” in the forms 


rp eint, logz = —if+ > 0g eing 
n=0 (3.1) 


where constant terms, which will not be needed, are ignored. Some 
of the early coefficients in these expansions are 
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and further coefficients can be found if desired. The expansions 
(3.1) may now be substituted in (2.3) and, assuming that orders of 
summations may be interchanged, the result can be written 


w= Agl+ s C,, emit > D, e~ ing, (3.3) 
n=1 n=1 
where 


n oo 
Cr at COn4y+0°B, Agt >-B, A,+ > "Ve B, (n => 1) 
r=1 r=1 
D, = cayt+ > 'y4B 
r=1 


D, = Sy» B, (n > 2) 


r=n 
The aerofoil 7 = 0 is now a stream-line if 
C,=D, (n>1). (3.5) 
The next step in the solution is to find the coefficients A,, and 
inspection of the equations indicates that it is probably possible to 
express A, as a power series in 1/b. We therefore assume that 
A, => %4,/* (n>9); (3.6) 


: k=0 
then, from (2.4), 


= 2,[a0%dat 2 ra0'Aajbr] [ores (3.7) 


and, assuming that the order of summations in the resulting double 
and triple series in (3.4) may be changed so that the series may be 
summed ‘diagonally’, the constants C,,, D, now take the forms 


C,, ce Gn4+iT DAL i°B,, °Ay+ by "By, "A; | /ok+ 


s "Yn m9 9A ml /b¥+- 
+2 1m= 


o k-1 ¢£ es 
~ > > “Yn Ope k-tA 
k=2 r=1 m=1 


x k 
Cag+ > Z m™y_4 i coat 
k=1 m=1 
oo r = 
$ SOS Sy may tA, mld, (8.9) 


k=2 r=1 m=1 


oo k 
0 k 
D, = Zz mY —n Kg °A —m/ > 
k=n m=n 
Ce) 


bk (n>1), (3.8) 


ory 


ae > May Mary, k-14 _.,/b* (n > 2). (3.10) 


k=n+1 r=n m=n 
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The boundary conditions (3.5) must be true for every value of b, 
so that, using (3.8)-(3.10) and equating coefficients of corresponding 
powers of 1/b, we obtain the relations 


ca,+1 °B, °A, +18, 14, = cay 
CO 44 +08, “Agt > "B, "Ag = 0 (n ae 2) 
r=] 


; (3.11) 


1Y_1 1g (n= 1), (3.12) 


n 
d "B, Mit 2 Bn i+ "Yn tag Ag 0 (n > 1) 


n 
iB, "4+ + 2. "B,, "A,+ $ "Yn "Og "Ay- m+ 


> “vo n Xo °A —m 3 
m=n = 
<= | (k>2;k>n), 
‘74. "dig'Ag (n = k > 2), 
0 (n>k>2) J 

From these equations the coefficients "A, for n > 1 can be ex- 
pressed successively in terms of the circulation coefficients °4,. Six 
of these coefficients are given below, but, apart from the length of 
the expressions, further coefficients could readily be obtained. Thus 





14 5 = CA9(Gy—a)—ia, “Ay 
‘A, $09(@2—AGy)"A y—ia,°A, 
"Ay = fed Gg(A9Gy— 2a, a+, G,)—ia, A,— 
— }aio(Gp—Ag)°A, + Hiag{ay Ag+ (@y—Ay)(4, —4,)}Ay 
= CAg{A, (Fy—Ay)—Ag Ag} —1(Ay y+ $07) y 
$4g{Ay 43 —Ay (Zy—A1q)}°A p—1(4g 4+ 3.03), 
= CAp{ (Gy—Ag)(Ay 4q+-A7) — 29 a, Ag—AG 4} — 
—1(a§ 43+ 2ay a, 4+ faz) A, 
With the help of (3.6) and (2.4) the expression (2.8) for the forces 
on the cylinder ey X = 0 and 


Yy 2) 
Se. _yrt14 J 
= cA, 20 > | ay syraigrr A +( ) rt 


S . (n+r)!i"- —f A, 4, 
i >< n—1)\(r— 1)! Q"+" prensa? (3.15) 
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so that, by writing © 7, = > Y,./b*, (3.16) 
k=0 


the coefficients Y; can be found once and for all for any given type 
of cylinder. The first few coefficients are 

Y,/p = 2c Ag, 

|p = 2ncA,—n7(%,)?, 
Y,/p = 2nc °A,—2n 4A, + 42c{a9(4)>—a_)+49(ay)—G,)}° 
—thin(a,—d,)(Ay)?, (3.19) 
Y3/p = —}1c*dyAo(Gy—4y)(4y—Gy)+ 
+ dric{(a, —G,)(2a9 Gy—y A, — Gy Gy) —02 A, +-2 G,}°A 9+ 

+-41c(2ay Gyp—Ay Ag—Ay 4y)°A, + 


| 


{a2 + G2 — 2a, dy — 29 ig +3(ay Ag+-Gy Gy)}(Ay)?— 


a 
—in(a,—d,)A_"A,—274,A,—77(%,)?+ 27%. (3.20) 
The expression (2.9) for the couple can be dealt with in a similar 


way. 


4. To complete the solution of the aerofoil problem it is necessary 
to obtain a value for the circulation. Adopting the usual procedure 
for aerofoil sections possessing a sharp trailing edge we choose the 
circulation so that the stream leaves the trailing edge smoothly. 
The aerofoil transformation (2.1) can always be chosen so that the 
trailing edge is the point £ = —v7 and the condition for finite velocity 
at this point is that dw/df = 0 for £ = —z. From (3.3) this gives 

a 
Ayt+id (—)"n(D,,—D,) = 0 (4.1) 
if the condition (3.5) is also used. If the expressions (3.6), (3.9), (3.10) 
for A, and D, are substituted in (4.1), the resulting equation must 
be true for all values of b, so that, equating to zero the coefficients 
of each power of 1/b, we obtain 
"Ay = t¢(@)—A) 
Ay = 0(*7-1 *Xo—*y-1 *%9)"Ao 
4.= is s (—)"*4n(™7_,, ™iXg 9A pm — Vn Xp "Ap m) + 
n=1 m=n 
k-1k-1 fr a 
+1 > 3 (—}* Hig (5 _., "Cin “y r ei Oe ™ Oh -r (4 


n=1lr=n m=n 
(k > 2) 
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From these equations and (3.14) the circulation coefficients °A4,. 
can be obtained successively and explicitly in terms of the constants 
a,, which define the aerofoil section. Results are given here for the 
coefficients “45, °4,, °A,. Expressions for further coefficients are 
necessarily more lengthy, but there is no theoretical difficulty in 
finding them, and, for a flat plate in which all the coefficients a,, are 
zero except a, and dy, the work is relatively simple. Since results for the 
flat plate have been obtained as far as terms in 1/b* by Tomotika and 
Havelock, there appears to be no point in reproducing these results 
here. It should be noticed, however, that the method used here, 
although having some connexion with Havelock’s work, appears to 
be somewhat simpler in concept. 

The coefficients 4, °A,, °A, are 


"Ay = ic(dy—4p), A, = 3¢(G)—ay)” 


|, (4.3) 
A, = fic{agGo(4,.—ag) + (@y—Ap)*(A, —A,) + 3g Ao(%)—A)} 


and, from (3.17)-(3.19), the coefficients Y), Y,, Y, in the expression 
(3.16) for the lift are 


Y) = 2aipc*(ay—ap), Y,/Yo = t(ag—Gp) 


= = $(Gy—A)(4, —G,) + fg Gy— }(@y—ay)?+ 
0 
Ay Ao(Gg—Ag) 


= }(2d9 ip — Gy Ay —Gy Gs) 
4(Gy—ao) 


Applications of these general results are considered in the next 


section. 


5. Fig. 1 shows a circular-are aerofoil ABC near a plane wall. 





Fie. 1. 


The mid-point of the are is B, and D is the mid-point of the chord 
AC. Since the theory developed above requires that the origin of 
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coordinates should be ‘inside’ the aerofoil, B is a suitable origin and 
the required transformation for the circular arc is 

+o -ee 6 be 

z+ia(1—e-2i%)et? ~ \y__ gif-2ia 


so that the coefficients in (2.1) for this problem are 


9 


a) = ae*9+®) sin x 
a, = —iae sin?x 

e * . . e ? 
— (—i)#e-&"—Dagid sin”-1y cos2a (n > 1) 


l = 2asin 2a 


a 


n 


where @ is the inclination of the chord AC to the wall, 4a is the angle 
subtended by the are AC at the centre of its circle, / is the chord 
AC. The height 6 of B above the ground is related to the height H 

of D above the ground by the equation 
b = H+4iltanacos8#. (5.3) 

From (5.2) and (4.4) we find 
Y, = zplc* sin(@+«)/cos « 
Y, = —4#l¥)sin(6+«)/cos « 
__ sin@sin(@+a) , sin(@—a) 5 cos 20 
BY, 8 COS a 64sin(@+«)  64cos*a 32 

For a flat plate Tomotika (10) found that, as the distance of the 
plate from the ground decreases, the value of Y/Y, always decreases 
when @ is large, but, for small @, Y/Y, first decreases and then 
increases. Calculations by the present writer for a circular-arc aero- 
foil (1) seemed to show that Y/Y, always decreases even for small 
values of 6. The reason for this can be seen by studying equations 
(5.4). The term Y,/lY, for « = 0 and small @ is small, but, as « in- 
creases, the effect of this term becomes more marked and would 
therefore tend to keep Y/Y, less than 1 for a greater range of values 
of 1/b. This is best illustrated by taking a numerical example corre- 
sponding to the work carried out previously. Thus, if « = 11° and 
G= &, 


(5.4) 





Y l 
— = 1—0-1404| - 
Y, ( 


and for a flat plate (a = 0) the result is 
Y eee wee 
— — 1—0-04358(—)+0-06392(—] +... 
} b b 


0 


)+o-osrsa (7) +. (5.é 
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6. The effect of thickness on the results can also be seen by con- 
sidering the symmetrical aerofoil of small thickness, which corre- 


sponds to 
ay = fU(1+e)e”, a, = plee®, a, = W(1—e)e”, a, = —}lee® 6.1) 
a,=0 (n>4) sal 


where ¢ is small and real and / is the chord of the aerofoil. For this 
case 


Y, = mplc?(1+-e)sin 6 
Y, 





*1 _ _4(14<)sin@ 
sto (6.2) 
at => as(2-+ 5e+ 3c?) + $(1+-€)? sin? 

0 


For small 6 the term Y,/lY, is slightly increased numerically by the 
small thickness, but the effect on the term Y,//*Y, is somewhat greater, 
so that we might expect that the thickness would accentuate the 
tendency for Y/Y, to decrease and then increase as //b increases. 
It therefore appears that, for small values of 6, the effects of camber 
and thickness work mainly in opposite directions so that the flat 
plate might be expected to give a reasonable approximation to actual 


cases. 


7. In this section the first steps are given for the solution of the 
problem of the general aerofoil (2.1) in a uniform stream ¢ which is 


bounded by the plane walls y = a, y = —b, together with a circula- 
tion K = 2k around the aerofoil. If 
d=a-+b, a = m(b—a)/2d, (7.1) 


a notation used by Havelock (2), then a suitable complex potential 
function w which makes the boundary walls stream-lines is 


w = cz+ik fogsinh a4 —logeosh & 4. ia] " 
20 (—)-! ds : 2 a as = Re 
ss 2, rae logsinh 5; + A, logeosh sit ia}. (7.2) 


Similar functions to those in (7.2) were used by Howland and 
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McMullen (4) for potential problems relating to groups of circular 
cylinders. Using the expansions 


(<1 @ 


(s—1)! dz 





where constant terms which are not needed are omitted, and where 


2n 

G5, [1 

ae +1" 2n 2n+1 pe 

ig = (—)* ? mttag = 0, 
nN \% 


"a. = (=p ‘et ° 


nB om 
s 


B = 16 = ttana, 


ib f d)2"-1 
fol) = (—fa +95) 





fon iaB) = i(—{1+84) "9 
2n+1 | | ds 


the complex potential (7.2) may be expressed in the form 
w = cz+iklogz+ y A,z2-"+ 5 Bs. 
n=1 n=1 


- —"Bo) "os Agt+"B, Ag 
2, — Sine, 5 


where qui 
s=1 

From this point onwards the solution follows that given in §2 
except that results are now expressed as a power series in 1/d instead 
of 1/6. For a flat plate, terms as far as those given by Tomotika and 
Havelock can be obtained without much difficulty, and, except for 
the length of the expressions, there is no theoretical difficulty in 
carrying out the calculations for the general aerofoil (2.1). 
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The infinite cascade of similar aerofoils in a uniform stream may 
also be handled by a similar method, starting from a complex 
potential function derived from the elementary functions used by 
Howland (3) for potential problems relating to an infinite row of 
equal circular cylinders. 


REFERENCES 


. A. E. Green, Proc. London Math. Soc. (2), 46 (1940), 19. 
. T. H. Havelock, Proc. Roy. Soc. A, 166 (1938), 178. 
. R. C. J. Howland, Proc. Cambridge Phil. Soc. 30 (1934) 315. 
. R. C. J. Howland and B. W. McMullen, ibid. 32 (1936), 402. 
. L. Rosenhead, Proc. Roy. Soc. A, 132 (1931), 127 
. 8S. Tomotika, Rep. Aero. Res. Inst. Tokio 7 (1933), 357. 
. —— ibid. 8 (1934), 115. 
. — ibid. 157. 
\ ibid. 13 (1938), 472. 
. 8S. Tomotika, T. Nagamiya, and Y. Takenouti, ibid. 8 (1933), 1. 
. 8S. Tomotika and I. Imai, ibid. 12 (1937), 421. 
ibid. 472. 
. W. F. Durand, Aerodynamic Theory (Berlin), 2 (1935), 91. 








NOTE ON AN ASYMPTOTIC FORMULA 
CONNECTED WITH 7-FREE INTEGERS 
By L. MIRSKY (Sheffield) 
[Received 24 January 1947] 
LET r, s be positive integers, and /,,...,4, distinct positive integers. 
We shall assume throughout that r > 2, and we define yu,(n) as 1 or 0 
according as 1 is or is not r-free.* 
One of the problems which I considered in a recent papert was the 
determination of an asymptotic formula for the sum 
F(z) = > p,(n+hy)...u,(n+k,). 
0<N<x 

I showed, in particular, that the error term in this asymptotic 
formula is of the form O(a‘/*+s-)+«), The object of the present note 
is to sharpen this result. 


Define the symbol A Ms 


as 1 or 0 according as the system 


8 
of congruences in ” 


n+k,=0(modn,) (l<v<s) (1) 
is or is not soluble. 
Furthermore, let {n,,...,n,} denote the lowest common multiple of 
By gooey Mge 
I shall prove the following 
THEOREM. As x>00, 
F(z) ae Ax+ O(a2ir+D+), 
where the O-constant may depend upon parameters other than x, and 
A= > Hep -wled gfe 
Pie A Oe Uo 


the infinite series being absolutely convergent. 


8. 


Remarks. (1) This theorem is only of interest for s > 2, since for 
s = 1 it reduces, in effect, to 
x 


trl) = t( r) 


* An integer is called r-free if it is not divisible by the rth power of any 
prime. 

+ ‘Arithmetical pattern problems relating to divisibility by rth powers’, in 
course of publication in Proc. London Math. Soc. This paper will be referred 
to as A.P. 


4. O(arr+) +e), 
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However, the more precise result 
x 
f(r) 


is well known and can be proved in a few lines. 
(2) It can be proved (A.P., Lemma 6, case g = 1) that 


D( pP \Keqy---s kg 
A = []p--2 yf |, 


2 Helm) = 5 + O(x"") 


Pp 


where D(p"|k,,...,k,) is defined as the number of integers amongst 
1,2,...,9" which are congruent (mod p”) to at least one of k,,..., k,. 
This shows, in particular, that A > 0 provided that, for every prime p, 


D(p" |k,,.-+5 k,) < Pp’, 


i.e. the numbers f,,...,&, do not contain a complete system of residues 
(mod p’). 

On the other hand, if D(p"|k,,...,4,) = p” for some prime p, then 
it is easily seen that F(x) = 0. 

(3) We may note that the theorem above enables us to improve 
the error terms in Theorems 2, 3, and 4 of A.P. to O(x/+)+), 

Notation. Here x, «, 8 are positive numbers; « is an arbitrarily 
small positive number. All other small letters denote positive integers. 
The O-notation refers to the passage x oo, and the O-constants 
depend at most upon e,7, 8, ky,...,k,, «,B; 
n,) denotes the highest common factor of 7,.,...,%5; 

apd 


denotes the number of systems of numbers 2, J,,...,/, 


a 


* <2, 
n+k,=n,l, (l<ov<s); 


. k,;x) denotes the number of systems of numbers 
N, Ay,..+, Ag, 0,,...,6, such that 


%<@, 
n+k,=ahb, (l<v<s), 
a. 
Q,...a, > x"; 


d(n) is the number of divisors of n. 
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Lemma 1. The system of congruences (1) is soluble if and only if 
(n,;,n;)\(k;—k;) (lL<ti<j <8). 
In the case of solubility the solutions form precisely one residue class 


(mod {n, 


This result is well known.* 


LEMMA 3. 


where K depends at most upon 8, k,,..., k,.- 


Proof. We first observe that the assertion is trivial when the 
system (1) is insoluble. Assume, then, that (1) is soluble. Hence, by 


Lemma 1, (n;.0;) < |kyp—kj| <i <j<e). (2) 


Now 1...”,/{ny,...,%,} can be written in the form of a fraction whose 
numerator and denominator are products of the highest common 
factors of two or more n,’s at a time. Hence, by (2), 


Ny...Ng 


{95,0005 Rs 
and the lemma follows. 


Lemma 4. L(x; ky,...,k,3«) = O(a--D+*) + O(gtir+d+e), 


Proof. The proof is by induction. 
Lie;k;c)= Fl = FF 2 


nx a, <(x+ ky)!” nsx 
n+k,=ajb, a,>2* n+ k, = 0(mod aj) 
a,>x* 


ial >» (= +00) = O(a-"-D)4 O(a), 


ay, 


a, <(a+h,)"!” 
a,>x* 


Thus the theorem holds for s = 1. Next assume that it holds for 


* See, e.g., A. Scholz, Hinfiihrung in die Zahlentheorie (Géschen), Satz 31. 
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some s > 1. Let 8 be some number whose value is to be fixed later. 
Writing a,...a,,, = a we have 

TAB} Beg 5025 Higgs tt) 
=A F1 pO F1 )=L+Lysay. 


nor 
n+k,=a,b, 
n+ ke+1=a541bs+1 


a>xz* 
a 


a” 
r r 
> 1 Y a | 
pine ee x*<a<rBls+ )/s Lseceg Moth 
pee, be | 


x*<a<rBl(s +1) 


Tr Tr 
hl “ts 
O Pam toe Rest) 


+1}>, by Lemma 2, 
{j++ 541} 


2* <a<rhs+ls 


a’ 
Asi 


Oz a iat + O(aPe+Dis+*), by Lemma 3, 


ae" 
gi—or—D+e) 4 O(aPis+Dis+e), 


> 1 ) 4 d(n+k,,1)} 


541 bs+1.=n+ kets 


Qy...dg> 28 
= Ofat L(x; ky,..., k,; B)} 
= O(xi-Br-D+2¢) + O(a2lr++2), by assumption. 
We therefore have 
TAB} By 50025 Han) * 
a O(al-ar-D+€) 4 O(g2lr+D+2¢) 1 O(al-Bir—-W)+2€) 4 C(a-Bis+iis+e)_ 


Ss ‘ 
Putting B — ——— we obtain 
&B rs+1 


L(x; sity a a) = O(ai-or-D+€) 4 O(g2lr+D+2e), 
and this proves the lemma. 
Proof of the theorem. We begin by noting that the series defining A, 


amely 
namely 12(4,)...u(a, 2 Aa 
, az} , 


is absolutely sadiadnea re Lemma 3. 
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Next, since p,(n) = > p(a), we have 
a\n 


Fe) = YT pla)..ea,) = SY w(u)--.n(a,) 


N<x aln+ ky 


a In +k, 


yoo Ag Sr" 


= ala) la) i 


Q,...ds<z"l" 


L) rors, by Lemma 4, 


"ss 


r r 
=2 }-(A,)..-4(4s) sees *) + Ofer), by Lemma 2. 
Q...ds<zilr {aj,-. , as} og, , 


1 


= Ar+ Oz 


atts by Lemma 3, 
Q...as>axilr 1:++s, 
= Az+ O(2' 2ir+I)+e)., 
In conclusion I should like to thank Dr. R. Rado for reading the 
manuscript and making a number of most helpful suggestions. 





A NOTE ON REAL QUADRATIC FORMS 
By J. A. TODD (Cambridge) 
[Received 12 February 1947] 


1. Ir A and B are two symmetric nxn square matrices such that 
B is non-singular, and if the roots of the equation |A—AB| = 0 are 
all distinct, then the quadratic forms x’ Ax and x’ Bx can be reduced 
by the same non-singular linear transformation x = Ty to forms 
y'A,y, y’B,y, where B, is the unit matrix and A, is a diagonal 
matrix whose elements are the roots of |A—AB| = 0. If these roots 
are A,,Ag,...,A,, then the matrix 7’ has for its columns a set of (properly 
normalized) vectors © such that (A—A; B)&™ = 0 (¢ = 1, 2,...,n), 
and the proof of the theorem turns on the fact that the n vectors £® 
are linearly independent.* If the roots of |A—AB| = 0 are not all 
distinct, such a reduction may or may not be possible. In order that 
the forms «’Az, x’ Bx should be capable of simultaneous reduction 
to diagonal form, it is necessary (and sufficient) that, corresponding 
to any p-fold root A; of |A—AB| = 0, it should be possible to 
find p linearly independent vectors £ which satisfy (A—A,; B)é = 0, 
or, in other words, that the rank of the matrix A—A,;B should 
be n—p. 

A familiar theorem states that, if A and B are real and if the form 
«’ Bx is positive definite, then a simultaneous reduction to diagonal 
form is possible. An essentially equivalent theorem is that any real 
symmetric matrix A can be reduced to diagonal form by a real 
orthogonal transformation. From this it follows that, if A; is a p-fold 
root of |A—AI| = 0, then the rank of A—A, I is n—p. It is difficult 
to find, in the literature, a simple direct proof of this result, though 


* This argument is, I believe, familiar, but a precise reference is not readily 
accessible. The proof of the theorem given, e.g. in Bécher, Introduction to 
Higher Algebra (New York, 1924), 167, is inductive. A proof on the lines 
mentioned above is given, essentially, by Lamb, Higher Mechanics (Cambridge, 
1920), 220. The very similar problem of reducing the equation of a 
collineation to diagonal form, in the case of distinct characteristic roots, is 
worked out on the lines suggested above in Aitken, Determinants and Matrices 
(Edinburgh, 1939), 89, from which the orthogonal reduction of a symmetric 
matrix is deduced—a problem essentially equivalent to the one considered 
above. 
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many references to it can be found.* It is an immediate consequence 
of the theorem just mentioned, and, conversely, implies this theorem. 
But the proof of the reduction to diagonal form of a real symmetric 
matrix by orthogonal transformation given in most elementary texts 
is inductive in character,} and this proof is open to the aesthetic 
objection that it does not correspond to the practical method of 
obtaining the reduction of a given symmetric matrix, which is based 
on the determination of a set of orthogonal eigen-vectors. The pur- 
pose of this note is to provide a simple direct proof of the fact that 
the rank of A—A;J is n—p where 2. is a p-fold root of the real 
symmetric matrix A. By making use of this result we can then 
obtain (in an obvious way) an alternative, non-inductive, proof that 
A can be reduced to a diagonal form by means of an orthogonal 
transformation, on the same lines as that which is available when 
the characteristic roots of A are all distinct. 


2. The proof depends on the following simple lemma. 

Lemma. If C is a matrix of rank r, then its r-th compound C® is 
of rank one. 

For, since C is of rank r, there exist non-singular matrices P, Q 
such that PCQ = D, 


where D is the matrix in which the first r elements in the principal 
diagonal are unity and every other element is zero. Hence, by the 
Binet-Cauchy theorem,{ 

POCO — Dm, 
But the leading element in D” is unity, while every other element 
is zero. Hence D” is of rank one, and so, since P” and Q” are 
non-singular, C” is also of rank one. 


3. We can now prove our main theorem. 


THEOREM. I[f A is a real symmetric n x n matrix, and if r, is a p-fold 
root of the equation ¢(A) = |A—AI| = 0, then the matrix A, = A—A, I 
is of rank n—p. 


* See, e.g., Levi-Civita, The Absolute Differential Calculus (London, 1929), 
205-7 ; Hilton, Homogeneous Linear Substitutions (Oxford, 1914), 58 ; Muir, The 
Theory of Determinants in the Historical Order of Development, vol. iii (London, 
1920), 294-5. 

t+ See Bocher, loc. cit., or Ferrar, Algebra (Oxford, 1941), 151-3. 

t See, e.g., Aitken, op. cit. 93-4. ; 
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If 6(A) denotes the kth derivative of 4(A) with respect to A, the 

conditions that A, be a p-fold root of ¢(A) = 0 are 

$A) = 90, g*A)=9O0 (L<k<p), gPA,) 40. (1) 
By successive application of the rule for differentiating a determinant 
it is seen that ¢™(A) is the sum of all the principal (n—k)-rowed 
minors of A—AJ, multiplied by the numerical constant (—1)*k!. 
Hence, using the notation of compound matrices, conditions (1) can 
be written in the form 

trace[Af!’-”] = 0 (0<k <p), trace| A{"-?)] + 0, 

since AY” is the 1 x 1 matrix whose element is |A,}. 

From (2) it follows that A{'-” is not the zero matrix. Hence the 
rank r of A, is n—p at least. Suppose, if possible, that r > n—p. 
Denote the matrix Aj? by C. Then C is a symmetric matrix (since 
A, is), and is real (since every characteristic root of A is real). If N 
is the number of rows and columns in C, then, since trace C = 0, 

Cy, t+Coot...+Cyy = 0. (3) 
By the lemma, C is of rank one. Hence, since C is symmetric, 
Ce = CS na Crs Cy = Cr Cos» 


and so, from (3), 
N N 


0= (Cy +Cyoot --. + Cyn)? _ > y Ore 


r=1 8= 


Thus, since C is real, every element of C vanishes, which gives a 
contradiction. Thus r = n—p and the theorem is established. 














THE SIMULTANEOUS REDUCTION OF TWO REAL 
QUADRATIC FORMS 


By W. L. FERRAR (Ozford) 
[Received 30 March 1947] 


1. Introduction 
Two real quadratic forms A(xz,2) and C(x,x), of which the latter 
is positive-definite, can be expressed by means of a real non-singular 
transformation in the forms 

> A,X? and > X?. 
This is commonly proved without finding the actual transformation 
from x to X. Such a transformation is not difficult to obtain when 
no two A’s are the same, particularly when C(z, x) is itself merely > 22. 
But there appears to be no readily available account of the actual 
transformation when C(x, x) is the general positive-definite form and 
the characteristic equation |A—AC| = 0 has multiple roots. 

The present paper gives the actual transformation under these 
general conditions. It is, essentially, an amplification and adaptation 
of the work of Levi-Civita in his Absolute Differential Calculus.* 
Levi-Civita’s work is not primarily concerned with a complete treat- 
ment of this particular problem, which is only incidental to his main 
purpose; it leaves out much that is essential to a full treatment of 
the subject. 

I have taken the occasion to use (§ 4) the technique of matrix 
algebra in place of the more usual methods of carrying out the 
transformations; once the algebra of the problem has been stated in 
terms of matrix equations, the actual reduction of A(x, 2) and C(x, x) 
to their canonical forms requires practically no calculation. 


2. The rank of the determinant |A—AC| 
2.1. In the field of real numbers let A, C be symmetrical matrices 
f ord . ae 1 - 
ot order 7 A= [a,;], o= [cis]; 
let A(x,x), C(x,x) be the corresponding quadratic forms in the n 
; and let C(x,x) be positive-definite. Then the 


n? 


equation |A—AC| = 0 


variables 2,, %p,...,X 


* pp. 205-7. 
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has n roots A,,Ag,...,A, (necessarily real) and there is a (real) non- 


mi fate 


singular transformation x = 7'X such that* 
A(x,x) = A, X?+A, X3+...+A, X2, 
C(x, x) = X3?+ X3+...4 X2. 


By suitable nomenclature we may take any one root of |A—AC| = 0 
to be A,. Let A, be an r-ple root. Then,} for any A, 


\a,;—Ac,;| = |T|?. L(A), 
where L(A) is the determinant of a diagonal matrix having A,—A in 
the first r places and A,,,,—A,...,A,,—A in the remaining n—r places. 
When A = A,, L(A) = L(A,), which has precisely r zeros in the 
diagonal and so is of rank n—r. Since T' is non-singular{ 


is also of rank n—r. Hence§ the equations 


> @5—Arcy&; = 9 (7 = 1,2,..., 0) 


i=1 
have r linearly independent solutions. This establishes the funda- 
mental result on which all the rest of the work depends; when A is an 
r-ple root of |A—AC| = 0, the equations 


n 


a (a,;—Ac,;)E; =0 (j oe 3. %....,%) 


i=1 


have r linearly independent solutions in the &;. 


2.2. When C(a,x) = > 2?, this fundamental result can be estab- 
lished by the direct and elegant method discovered by J. A. Todd.|| 
It can then be extended to the general positive-definite form C(x, x) 
by using invariance arguments similar to those given in § 2.1. 


3. The ‘principal directions’ of two quadratic forms 
3.1. We use the summation convention** and write our quadratic 


forms as A(x, x) = a,;x'x', C(x, x) = ¢,,; x*2/, 

* Ferrar, Algebra (Oxford, 1941), 151. Other references to this book will 
be indicated by F., followed by the page number. 

t F. 127, 8. t F. 110. § F. 101. 

See the preceding article. The result of § 2.1 is not new, but it is not easy to 

find, in the literature of the subject, a satisfactory proof that is free of 
difficult ‘invariant-factor’ arguments. 

*e F. 4i. 











188 W. L. FERRAR 


the n variables being a* (i = 1, 2,...,). We-use the single letter 2, 
when convenient, to denote the x variables 2‘. We take C(x,x) to 
be positive-definite. 
We shall find n sets of values 
=i. (r = 1,§,...:0) 
for which the ratio 
a,; X}+-a,, X?+...+a,;X? 


+@_.. 
2 Ms. 3 


1 a 72 yn 
Cy; X}+-Ce;, X2+ ... 6,5; XF 


is independent of j, or, on using geometrical language, the vector 
with components 
Gm Bin MS, «sci Mites 
is parallel to the vector with components 
a Aa oe (me oe Ae 
For convenience, we call the latter a ‘principal direction’ vector of 


the two forms; the vector 


1 re é 
x1, x2, ..., X” 


r 


from which it is derived we call a ‘A-vector’. 


29 Wass, ’ - . —f. 
3.2. Two different roots; A, ~ , 
Let A, be any root of the characteristic equation 
Aj;—ACj; = 0. 
Then there is a non-zero A-vector 
ra r 2 7n 
X}, 4LaA» eeey ah 


rar urht * 4) ss + 9 
for which (a,;—A, CX) = 9 (jf = 1,2,...,m), 
and, since only ratios of the X are involved in (2), the components 


can also be made to satisfy 
ot 
Cry Xiq) X4y) = 1. (3) 


It is at this point that the positive-definite character of C is essential 
to the argument. 

Let A;, ~ A;, be any other root of (1) and let X,, be a corresponding, 
non-zero, A-vector. Then 


(4;—Ay. C4) X ix) = @, (4) 


* When a suffix is enclosed in brackets no summation with respect to that 
suffix is implied, even though the suffix is repeated in the formula. 
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Multiply this by X} and (2) by Xj, sum for j = 1, 2,...,n, and subtract 

the two results; we get, since the matrices A and C are symmetrical, 
(A, —Ax)(Ciy Xj, X}) = 0. 

Since A, + A,,* c,; Xj, Xi, = 0. (5) 


3.3. A multiple root 
Let A, be a root of (1) of multiplicity r. Then, by § 2, there are 
r linearly independent solutions of (2), each of which can be made 
to satisfy (3); and each of them also satisfies (5) for any A, + A,,. 
When A,, is a double root of (1), let 
Xi,, Xi, = 1,3....,2) 
be any two independent solutions of the equations (2), each satisfying 
(3). Then 
either Ci; Xian Xduje = 0, 
or we can replace X,. by X,.+6X,,, choosing 6 so that 
C45 Xia (Xbyet PXfay) = 0; 
for in the last equation the coefficient of @ is unity, by (3). We then 
re-name the solution X,.+6X,,, calling it, for convenience, Xj». 
With this nomenclature, 


Cj Xian Xlnpe = 0. (6) 
When A, is a triple root of (1), there are three linearly independent 


sets of solutions of (2) and we proceed in the same way. Having 
chosen the first two sets of X to satisfy (6), either the third set X,, 


satisfies ; ; : : 
Ci; Xia Xinys = 9, Cij Xin2Xims = 9, 
or we can replace X,,, by 
XpgteXpetvXj1, (7) 

choosing » and v so that 

C45 Xtuya(Xtnyst+eXlnetrAtan) = 0 
and Ci; Xinyo( Xing + ue XGpyot+vXiq)1) = 0. 
In view of (3) and (6) these reduce to 


> 


ie i Xi ce i YF 
v= —C; Xan Xis B= —Cy Xie Xs 
Thus, using (if necessary) the notation X,, to represent what we have 


* The work of § 3.2 is, of course, not new; all the ideas are familiar and the 
work is included here so that the proof given may be complete in itself. 
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previously denoted by (7), we have a set of three linearly independent 
A-vectors X;,,, Xj2, Xj3 satisfying 
C4; Xia Xt) = 1, (8) 
C1; Xtwr Xins =0 (r #8). (9) 
Similarly, when the equation (1) has an m-ple root A,, we can find 
m sets satisfying (8) and (9). As we saw at the beginning of this 
sub-section, each of these also satisfies (5) for any A, + A,. 


3.4. Restatement of the results 
We may sum up in a more convenient form thus. Let A,, Ag,...,A,, 
be the » roots of (1). We have found n A-vectors X,, Xg,...,X,, 
satisfying the conditions 
C4; Xl) Xf) — 2. (10) 
Cy XX} — 0 (r af 8). (11) 
We may call such a set of A-vectors a ‘A-frame’ and each member 
of the set a ‘frame vector’. To each simple root A of (1) there corre- 
sponds a unique frame-vector (§ 2.1, with r = 1). To each m-ple root 
there corresponds a set of m frame-vectors, which is not unique, but 
allows an infinity of choice. For example, with a triple root and any 
given three independent solutions of (2), say 
Xn Xz» Xs: 
we may select an arbitrary 
Xi t OX ot $Xis 
to be our initial X,, of § 3.3. The situation is the familiar one met 
with in determining principal axes for a conicoid of revolution. 
Finally, to each frame vector X,. corresponds one principal-direction 
vector &,, given by (ef. § 3.1) 
G=GyAt (6 > 250009 M). 12) 
4. The transformation 


4.1. We first express the results of § 3 in matrix notation. Indi- 
cating a matrix by the element in its rth row and ith column, let 


Az=[a,), C=[egl, 
X = [X], € = (€]. 


Then the equations (12) of § 3, defining € in terms of X, become 
= X'C, (1) 
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while the ‘orthogonal relations’, (10) and (11) of § 3, become 
X’CE = I. (2) 
From (1) and (2), €X = ff, €&= I-1, (3) 
Finally, the equations (2) of §3, satisfied by the A-vectors, become 
C’'XA = A’E, (4) 


where A is a diagonal matrix with elements A,, Ag,...,A,,. This results 
from the following considerations. 

When a square matrix is post-multiplied by a diagonal matrix 
with elements A,, Ag,...,A,, its ith column is multiplied by A;. Accor- 


dingly, when O'X = [ey X3] 
is post-multiplied by A, the product is* 
[Ata Cor X71], 
which, in virtue of equations (2) of § 3, is 
[a,,X9] = 4’. 
This proves (4). 
4.2. The transformation of the two forms 


The transformation that reduces A(z,x) and C(x,x) to their 
canonical forms is 


y,= 2 (ry = 1,2....,2), (5) 
or, in matrix notation, gy = &e, (6) 
where y and x are single-column matrices with elements ¥;, Y2,.--,Yp 
and 2,,%,...,%,. This transformation is non-singular, as (3) shows. 
Now y’y is a matrix having the single element 
Yityet...+y%- 
But, by (6) and (1), y = 2 = I'Cz, 
and y = 2'C's. 
Hence y'y = «'C' Xéz, 
or, on using (3), y'y = 2'C's. 
Again, y’Ay is a matrix having the single element 


Ay Yi AL YE tn Un 
Now, by (7) and (8), 


y’Ay = 2'C'X Aga = x'A'XEx, 


* The bracket about the suffix i in the formula indicates that there is no 
summation with respect to 7. 
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on using (4). But € = X-1, by (3), and so 
y’ Ay = 2'A’x. (10) 
The matrices on the right of (9) and (10) are the matrix expressions 
for the quadratic forms* C(x,x) and A(zx,2). Hence the transforma- 
tion y= Ex (r = 1,2....,%) 
yields the identities 
C(x, 2) = yit+yit+...t+yr 
A(x,x) = A, y? +A. y3+...+A, y2. 


* F. 122. The matrices A and C are symmetrical, so that A’ = A and 
C’ = C. We could, of course, use rows instead of columns in forming the 
equations at the end of § 2.1. 
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